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ABSTRACT 

This thesis presents an application of a self-tuning 
regulator to a pilot plant double effect evaporator. A digital 
computer simulation study was performed to investigate the effects of 
several design parameters and non-ideal operating conditions on the 
performance of the regulator. Experimental runs were then made using 
an IBM 1800 process control computer to evaluate the performance of the 
self-tuning regulator and to verify the simulation results. 

The digital simulation results revealed that the self-tuning 
regulator performed better than a conventional controller when process 
noise and unmeasured disturbances were considered. They also demon- 
strated that the regulator can generate a satisfactory set of control 
parameters even though poor initial estimates are used. The effects of 
using different design parameters in the self-tuning regulator were 
studied and the results demonstrate that the regulator gives satisfactory 
control even for a poor set of design parameters. It was also found 
that the initial values of the design parameters can often be improved 
after examining the output response, the control variable response, and 
the parameter estimates. 

In general, the experimental application was successful eyen 
though some difficulties were encountered in the experimental runs for 
conditions that were similar to those used in the simulation study. 
This was due to the interactions with other control loops in the 
evaporator system. The controller gains in the level control loops had 
to be decreased before the self-tuning regulator would function 
properly. To summarize, the performance characteristics of the self- 


tuning regulator were similar in both simulation and experimental] 
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studies and the successful application of the regulator to the 


evaporator system was demonstrated. 
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CHAPTER ONE 
INTRODUCTION 


Automatic control system design usually requires an accurate 
model of the system to be controlled before a suitable controller can 
be designed. However for many industrial processes, it is very 
difficult to get an accurate model because of changing process 
conditions and the disturbances in the system. Furthermore, even though 
identification schemes can be used to determine models for stationary 
processes, they may involve an excessive amount of plant experiments and 
off-line analysis. Thus an adaptive control system in which the 
parameters of a simple model are estimated on-line and subsequently used 
in the controller is very attractive for practical applications. 

In the self-tuning regulator [1], the parameters of a single- 
input, single-output stochastic model are estimated on-line using a 
linear least squares recursive technique; the updated parameter 
estimates are then used in a minimum variance controller to calculate 
the control signal at each sampling instant. Very little a priori 
information about the process is required and the regulator can be 
easily implemented. This approach takes into account system time delays 
and can be used to control processes with slowly time varying parameters. 
The self-tuning regulator, as developed by Astrom and Wittenmark [1] and 
Peterka [6], is the chief concern of this thesis. 

This thesis presents an application of the self-tuning 
regulator to two systems: a simple stochastic model with constant but 
unknown parameters and the pilot plant double effect evaporator in the 
Department of Chemical Engineering at the University of Alberta. The 


evaporator application includes an extensive digital computer simulation 
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plus experimental verification on the actual pilot plant. 


1.1 Objectives of the Study 
The chief objective of this thesis is to evaluate the utility 


of the self-tuning regulator for practical applications through 
simulation studies and experimental application to the pilot plant 
evaporator. The self-tuning regulator has been previously studied in 
several simulations and industrial applications. However, most of the 
studies considered only single-input, single-output systems. In this 
thesis, the self-tuning regulator is used in the composition control loop 
of the evaporator which is an interacting multivariable system. 

A second objective of this thesis is to investigate the effects 
of various design options on the performance of the self-tuning 
regulator, and also the effects of "realistic" process conditions. An 
on-line experimental application of the self-tuning regulator is desir- 


able to verify the simulation results. 


1.2 Structure of the Thesis 

The thesis consists of three major sections. The first 
section, Chapter Two, begins with a literature survey and a deriyation 
of the self-tuning algorithm. It also includes simulation results for 
a simple numerical example selected from the literature to illustrate 
the behaviour and characteristics of the self-tuning regulator. This 
Simulation example also provides some familiarity with the design 
parameters of the regulator. 

The second section, Chapters Three and Four, presents the 


simulation results for the evaporator system including the effects of 
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different choices of the design parameters that are used in the self- 
tuning regulator. These results were useful in developing practical 
guidelines for the choice of design parameters in the experimental study. 
This section also includes a description of the actual evaporator and 
the IBM 1800 process computer. 

The third section describes the experimental application of 
the self-tuning regulator to the evaporator system and provides a 
verification of the simulation results. This material is presented in 
Chapter Five. 

The general conclusions drawn from this investigation are 


summarized in Chapter Six. 
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CHAPTER TWO 


MINIMUM VARIANCE AND SELF-TUNING REGULATORS 


2.1 Introduction 

For a single-input single-output stochastic system with known 
constant parameters, the corresponding minimum variance regulator can 
be easily derived. This regulator will have constant parameters and 
will minimize the output variance. However, in most industrial processes 
that can be represented by mathematical models, the parameters of the 
models are either not accurately known or are time varying due to 
changing process conditions. Consequently, in recent years there has 
been considerable interest in control strategies in which the regulator 
parameters can be autcm>*icztty adjusted on-line to compensate for 
unknown process dynamics or changing operating conditions. Such a 
regulator is often referred to as an adaptive controller, and one 
class of such controllers that will be studied in detail in this 
thesis is the self-adjusting or self-tuning regulators proposed by 
Astrém and Wittenmark [1]. This chapter presents a literature survey 
and the basic equations for the minimum variance and self-tuning 
regulators. Simulation results for a numerical example due to Astrom 


and Wittenmark [1] are given at the end of the chapter. 


2.2 Literature Survey 


System identification techniques for different classes of 
mathematical models have been extensively studied in recent years. 


This interest is reflected in the survey papers by Astrém and Eykhoff [2], 
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Nieman et. al. [3] and Gustavsson [4]. From an industrial viewpoint, 
a chief objective of process identification is to obtain better 
knowledge about existing plants so that their operation and control 
can be improved. Therefore, in many applications the purpose of 
identification is to develop a mathematical model that describes the 
process well, and can be used to design a control system for that process. 

As noted in the previous section, there has been a large 
amount of interest in the development of adaptive control strategies 
which combine estimation and control and can be used on-line on a 
real-time, process control computer. The approach of interest in this 
thesis, self-tuning regulators, stems from stochastic optimal control 
theory. It can be shown that if the design of a control system is 
approached through identification then according to the "separation 
hypothesis" [2,8], the optimal regulator can be divided into two parts, 
identification of model parameters and their incorporation into the 
optimal control law. 

The stochastic transfer function models and the real-time 
estimation techniques used in the self-tuning control strategy were 
first discussed by Astrém and Bohlin [5]. Peterka [6] was one of the 
first researchers to develop the complete estimation and control 
algorithm. The most extensive theoretical, simulation and experimental 
work on the self-tuning regulator has been performed by Astrom and 


Wittenmark [1,8], Wittenmark [7], and Astrém Sbacal alo te 


de 


2.2.1 Basic Theory 


Early papers on stochastic identification and control 
methods were mostly concerned with separate identification and control 
steps. That is, generally an off-line identification scheme would 
be used based on experimental input/output data. The identification 
step would be followed by the design of a control strategy and 
subsequent implementation on the system. Identification problems are 
generally characterized by three items: the class of models to be 
used, the class of input signals, and a criterion for the estimation 
[2]. Therefore, for the identification phase, certain parametric 
models and the type of input signals that will generate consistent 
parameter estimates have to be selected. Also, the minimization of 
a scalar loss function is generally specified as the criterion for 
estimation. Once the identification phase is accomplished, the design 
of a suitable control strategy can be achieved using existing control 
theory. This type of separate estimation and control for a single 
input/output etoBeStic system has been studied by Astrém and 
Wittenmark [8]. 

The idea of combining the estimation and control steps was 
considered by Kalman [10] in 1958; subsequently, there has been a 
large amount of research on this topic. The type of combined estimation 
and control method that is considered in this thesis is the self- 
tuning regulator. The regulator is designed based on the assumption 
that the system can be represented by a single-input single-output 
discrete stochastic model of the class considered by Astrom and Bohlin [5], 


and which is generally referred to as the Astrom-Bohlin model. Peterka 
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[6] used this model and derived an algorithm which automatically 
adjusted the parameters of a digital controller on-line. The basic 
structure of the algorithm can be seen in figure 2.1. It consists 

of two parts: real-time recursive least square estimation of model 
parameters and a minimum variance control law for the estimated model. 
Peterka has shown that under certain conditions, the controller is 
the same as the minimum variance controller that could be derived if 
the system parameters were known. 

Astrém's involvement in the development of the self-tuning 
regulator began quite early. He and Bohlin [5] presented a general 
form of stochastic models and investigated the identification of linear 
Stochastic systems, using the maximum likelihood method. In a later 
paper [26], Astrom examined the achievable accuracy of the parameter 
estimates and showed that his results can be applied to generalized 
least squares and model adjusting techniques. Astrom and Wittenmark [8] 
also did a detailed study of separate and combined estimation and control; 
they used models with constant but unknown parameters and models in 
which the parameters are time varying Gauss-Markov processes. Their 
minimum variance control law uses previous input/output data, the 
estimates of the parameters and the variances of the estimation errors. 
This method of combining least squares estimation and minimum variance 
control law was simulated using models with unknown, time varying 
parameters by Wittenmark and Wieslander [11]. 

Astrom and Wittenmark [1,7] then applied this type of control 
strategy to models with constant but unknown parameters. They assumed 


that the system was minimum phase with known time delay and that a 
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Figure 2.1 Basic structure of self-tuning regulator . 
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bound can be given to the order of the system. They also proved two 
theorems for the closed-loop system that are valid if the parameter 
estimates converge. (Unfortunately, general conditions guaranteeing 
convergence have not been reported, though efforts were made in this 
area.) The first theorem states that under weak assumptions, the 
covariance of the output variable and the cross covariance of the 
control and the output variables are both zero. In the second theorem 
it is shown that the control law obtained converges to the minimum 
variance control law that could be derived if the constant model 
parameters were known. 

Similar self-tuning regulators have been designed for the 
same stochastic model using different estimation techniques such as 
stochastic approximation [12] and maximum likelihood [13]. Other 
regulators have been investigated by retaining the least squares 
estimation part but using a different performance index which minimizes 
the output variance and penalizes the control signal. Studies on this 
type of regulator has been reported by Cegrell and Hedqvist [14] and 
Clarke and Gawthrop [15,16]. 

The asymptotic properties of the self-tuning regulator 
have been examined by Ljung and Wittenmark [17], they have presented 
the techniques and basic theorems to analyse this problem. They have 
shown that for closed-loop systems near or outside the stability 
boundary, the regulator stabilizes the system even if the model noise 
does not agree with the true noise characteristics. They have also 
shown that the regulator does not converge for general noise structures, 


and they have constructed an example to illustrate this. 
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The text by Astrém [13] provides an excellent introduction 
to stochastic control theory for linear systems with quadratic per- 
formance criteria. Its analysis of stochastic systems, parameter 
estimation and optimal stochastic control includes both the continuous 


and discrete cases. 


2.2.2 Modification and Extensions 
Wittenmark [7] has extended the self-tuning algorithm 
to include feedforward control for known disturbances and setpoint 
control of the output variable. 

For non-minimum phase systems, the standard self-tuning 
regulator cannot be applied in general. This is because the minimum 
variance controller is extremely sensitive to variations in the 
parameters [1,13]. Several methods have been proposed to overcome 
this difficulty. Astrom 13d} Astrém and Wittenmark [18], Turtle 
and Phillipson [19] and Peterka [6] have suggested the use of sub- 
optimal control strategies. Clarke and Gawthrop [15,16] claim that 
this problem can be easily solved by including the square of the 
control signal in the quadratic performance index. 

An extension of the self-tuning algorithm to multi-input 
systems can be made quite easily by simply adding the appropriate terms 
to the Astrém-Bohlin model and the algorithm [23]. However, an 
extension to multi-output systems is much more difficult since a new 
general structure of stochastic models has to be formulated. Further- 
more, the complexity of the algorithm and the computational requirements 


will increase tremendously. Gustavsson, Ljung and Sdéderstrom [20] have 
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reported their theoretical analysis on the identification of linear 
multivariable process during closed-loop operation. Astrom and Peterka 
[21] have considered the extension of the self-tuning algorithm to 
multivariable systems with constant but unknown parameters. However, 
the extensive on-line computational requirements include an orthogonal 
transformation and the solution of a matrix Ricatti equation each time 
the parameters are updated. Consequently, their multivariable self- 


tuning regulator is not attractive for most practical applications. 


2.2.3 Industrial Applications 


Several applications of the self-tuning regulator to 
industrial processes in Sweden have been reported [9]. Wittenmark and 
Borisson [22] have used the regulator to control the moisture content 
of paper from a paper making machine. They reported that very good 
results were obtained particularly during the startup phase and for 
steady state control. They found that the number Peaite parameters 
in the regulator was not very crucial in their experiments since it 
had little influence on steady state performance. They also found 
that the regulator can handle non-stationary disturbances without 
difficulty. Feedforward compensation of a measured disturbance, couch 
vacuum, was included in the algorithm. 

Cegrell and Hedqvist [14,23] also used the self-tuning 
regulator to control a different paper machine. They considered the 
moisture control loop of the machine which consisted of two inputs; 
thick stock flow and steam pressure, and one output; the moisture 


content. They found that replacing a conventional digital PI controller 
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by the self-tuning regulator resulted in a decrease in output 
variance and a reduction of losses during quality changes and mil] 
Setup. They also confirmed that the assumed order of the process 
model is not important. 

The application of the self-tuning regulator on a mine ore 
crushing machine has been reported by Borisson and Syding [27]. For 
changing characteristics of the incoming ore and the crusher itself, 
the regulator was able to adapt to the varying situations and control 
the input ore flow to the crusher to maintain high production. Their 
results demonstrated a 10% increase in production in comparison with 


a fixed parameter PI controller. 


2.3 Minimum Variance Regulator 


Consider the Astrom-Bohlin model which is a single-input 


single-output stochastic model of the form [13]: 


y(t) + ayy(t-1) eae ay (t-n) = bju(t-k-1) + ee t b ,u(t-k-n) 


+ ALe(t) + c,e(t-1) eet c,e(t-n)] i251) 


where 
y is the scalar output 
u is the scalar control 
t is the sampling instant 
k is the time delay (a non-negative integer) 
n is the order of the system 


X» is a constant 


and {e(t)} is a sequence of independent, normal (0,1) random variables. 
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Equation (2.1) can also be expressed as: 


A(q7!) y(t) = B(q7!) u(t-k-1) + ac(q7!) e(t) (2.2) 
or 
-1 al 
y(t) = B(q__) u(t-k-1) + 2 C(q_) (t) (253) 
a(q"') eS to 


where qu! is the backward shift operator, i.e. qu! y(t) = y(t-1) and 


] 


A(q7!) =] + aq +... tag 


B(q7!) Saba boa" +... + bat 


C(q7!) 1 + ae +... + cq 


Before the minimum variance control law can be derived, certain 
conditions have to be imposed on polynomials B(q7!) and C(q7'). The 
polynomial B is normally assumed to have all its zeroes inside the unit 
circle, i.e. the system of equation (2.1) is minimum phase. Polynomial 
C is also assumed to have all its zeroes inside the unit circle. These 
conditions are explained in detail by Astrom [13]. 


Equation (2.3) can be written as 


=) - 
y(t+k#1) = Neh 3 u(t) +2 sai : e(t+k+1) (2.4) 
q q 


In order to derive a minimum variance regulator it is convenient to 


express C(q7!) as [1]: 


Games Alqabyeztael) Hig Sh atqnh) (2.5) 


Z(q7!) Rlehe2,G tes i 2g (2.6) 
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Introducing equation (2.5) into equation (2.4) gives 


-] hy a eee 
a(q") a(q') ; 
(2.8) 
or 
-1 -] 
y(ttket) = BUC) ycey + az(q7!) e(teket) + a SQ) e(t) (2.9) 
A(q") A(q ) 
But from equation (2.2) 
_ AC!) B(q7!) -k-1 
ne(t) = “AT! y(t) - (t) (2.10) 
C(q ) c(q7!) : : 


Replacing Xe(t) in earstion /? 9) gives 


al -k-1 -| -| 
y(ttk+1) = aZ(q7!) e(t#k#1) + Hane EIeeal re 
Bide iuiiengy ).Cig .) 


aE) 75 (2.11) 
c(qt) > . 


which can be simplified to, 


2 a G 
y(ttk+1) = aZ(q7!) e(texe1) + BUG 7 a Du(t) + SO) yey 
C(q C 
(?. 12) 


A minimum variance regulator is a regulator that minimizes 


the output variance (or loss function) V: 


V = ELy*(t)] (2.13) 
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where E denotes the expectation operator. Let u(t) be an arbitrary 
function of y(t), y(t-1) ... and u(t-1), u(t-2) ... Then from equation 
(2°12): 


ELy“(t+k+1)] = E[aZ(q7!) e(t+k+1)]? 


2 
=] “4 = 
E BL Quel int g 2.) t G(q_) 0 | 2.14 
Re wed Fete | ce 


Since e(ttl), ... , e(ttk), e(ttkt+1) are independent of y(t), y(t-1), 
. and u(t-1), u(t-2), ..., hence 


2 2 


2 
( ni rezaibter, + 0] F215) 


Ey" (t+k+1) > 201 +.Z 


where the equality holds for 
B(q7') Z(q7') u(t) + (q7') y(t) = 0 (2.16) 


Re-arranging equation (2.16) gives the desired minimum variance control 


law: 


a 
ee) = ge dy (2) (2.17) 


If the minimum variance regulator is used to control the 
system represented by equation (2.1), the output would be a moving 


average process of order k, i.e. 
y(t) = 2z(q7!) et) 


or 


y(t) = ALe(t) + z,e(t-1) ye gote z,e(t-k)] (2.18) 


; ae 


2.4 Self-Tuning Regulator 


If a dynamic system has constant but unknown parameters, 

One way to formulate a control law is to first estimate the model 
parameters using some identification technique, and then derive the 
appropriate control strategy off-line by optimizing a selected per- 
formance criterion. This approach consists of separate identification 
and control steps. 

Astrom and Wittenmark [1] have proposed a type of adaptive 
control strategy that combines the problems of identification and 
control. In their self-tuning regulator, the parameters of the assumed 
model are updated every sampling instant using current and previous 
input/output data. The control signal is then calculated from a 
control law which utilizes the updated model. One desirable feature 
of their self-tuning regulator is that if the parameter estimates 
converges, the regulator converges to the minimum variance control 
law that could be derived if the model parameters were known. 

In the Self-tuning regulator, a recursive least squares 
algorithm is used to update the parameter estimates at every sampling 
instant, and then the control variable is calculated from a minimum 
variance control law which uses the updated parameter estimates as if 
they were the true parameters of the system. Astrom and Wittenmark [1] 
have noted that the estimation of the model parameters in equation (2.2) 
can be made Simpler if it is assumed that c(q7!) = 1, since then the 
standard linear least squares analysis can be applied. Therefore, the 


Simplified system model is written as 


A(q7!) y(t) = B(q7!) u(t-k-1) + ae(t) (2.19) 
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To simplify the on-line computation of the control law, the identity 
in equation (2.5) is incorporated into this model. Substituting 


c(q7!) = ] into equation (2.12) gives the “predictive model": 
y(t+k+1) - a(q7!) y(t) = B(q~!) z(q7!) Tig Me ,Z(q7!) e(ttk+1) (2.20) 


Equation (2.20) can be expanded to 


y(ttk+1) + oy(t) +... + a y(t-n+1) = B,Lu(t) + B,u(t-1) 
Sr ig: 2 B u(t-2) J + ¢(ttk+1) (2.21) 


where 2 = ntk-1 and c(t) is a moving average of order k of the driving 
noise e(t). The coefficients Os and B. are related to the parameters 


a. and b. in equation (2.1) through identity (2.5) and control law (2.17). 
If the coefficients are constant and known, then the minimum variance 


controller for the predictive model (2.21) is [7]: 


u(t) = a. [a,y(t) ae ek we ay (t-n+1)] - B,u(t-1) -..e - B u(t-£) 
0 
(2.22) 


If the parameters {a} and {85} of the predictive model are 
estimated on-line and the estimates, fa, (t)} and {8.(t)}, are used in 
the feedback controller (2.22), then the estimates are correlated with 
the disturbances {e(t)}. As a result, it may not be possible to 
accurately determine all parameters. To avoid this possible difficulty, 
Astrom and Wittenmark assumed that parameter Bo is constant and given. 


Then equation (2.21) can be written as 
y(ttk+1) + a(t) y(t) +... +, y(t-n#1) =  [u(t) + 8)(t) u(t+l) 


+... + é , (t) u(t-2)] + ¢(ttk+1) (2.23) 
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Next, introduce the column vectors, w(t) and 6(t): 


ct 
— 
u 


[-y(t) -y(t-1) ... -y(t~-n+1) B u(t-1) hap B u(t-2)] 


al 


t) = [a,(t) a(t) ... a, (t) AY ee 8 (t)I 


where superscript T denotes a vector transpose and the underline 


denotes a vector. Then equation (2.23) can be written as 


y(ttk+l) = 60 u(t) + w'(t) O(t) + e(tek+) (2.24) 


The least squares recursive equations that are used include exponential 


weighting of past data. These equations are [7]: 
A(t) = (t-1) + K(t)y(t) - 8) u(t-k-1) - yl (t-k-1) 8(t-1)] (2.25) 
K(t) = P(t) plt-k-1) [R + yl (t-k-1) P(t) w(t-k-1)]° (2.26) 


peter) = 1 gpce) - x(t) ER + yT(tek-1) E(t) v(t-k-T)] KU(t)? (2.27) 


|m> 
Ft 
ctr 

— 
— 


is the estimated value of o(t) 


K(t) is the gain vector {(ntz) x 1} 


Imo 
— 
cr 
ww 
— 


is the covariance matrix {(ntz) x (n+z)} 


is the exponential forgetting factor 


5 = 
— 


R is the covariance of the noise variable, e(t) 


The derivation of the recursive form of the least squares 
estimation algorithm can be found in references [6], [24] and [25]. 
The theoretical reasons for using the exponential forgetting factor 


are twofold. First, the rate of change of the parameter estimates 
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depends on the gain vector K(t) in the recursive least squares 
equations, and after a period of time the gain vector elements may 
have decreased to small values before the estimates reach the true 
parameter values. Placing more weight on recent data and less on 
past data tends to prevent the elements of K(t) from becoming too 
small. Second, for processes with slowly time varying parameters, it 
is possible to follow these parameters by using a forgetting factor 
of less than one. 

To summarize, the self-tuning regulator estimates the unknown 
parameters, fa. } and {8.35 on-line and these estimates, {a.(t)} and 
{8.(t)}, are then used in the minimum variance controller (2.22) as 
if they were the true values. The following calculations are performed 
on-line at every sampling instant: 

(a) Parameter estimates ta. (t)} and {8.(t)} are calculated based on 
input/output data using the recursive least squares algorithm in 
equations (2.25) to (2.27). 


(b) The control signal is calculated from: 
u(t) = 4 [a4 (t) y(t) +... + &(t) y(t-nt1)] 


ghey it )pult-1).-.....- B u(t-2) (2.28) 
which can be written as 


u(t) = 2 8" (t) v(t) (2.29) 


° 
Astrom and Wittenmark [1] have proved that if the parameter 
estimates converges, the control law converges to the minimum variance 


control law that can be determined if the parameters of the system are 
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known. This occurs despite the fact that the parameter estimates are 
biased due to model simplification and the correlation between parameter 
estimates and disturbances, 

In general, the least squares estimation technique wil] 
give unbiased parameter estimates under the following conditions [1,2]: 
the predictive model has the correct form; the residuals {e(t)} are 
independent; the input is persistently exciting; and the input 
sequence {u(t)} is independent of the disturbance sequence {e(t)}. 
However in the self-tuning regulator several of these conditions are 
not satisfied since the predictive model in equation (2.20) does not 
have the correct form and the inputs and disturbances are correlated 
due to feedback control. Furthermore, the input, u(t), is not necessarily 
"nersistently exciting". For these reasons, the parameter estimates 


generated by the self-tuning regulator tend to be biased. 


2.5 Numerical Example 


Before one can use the self-tuning algorithm, several constants 
have to be specified for both the estimation and control parts. The 
constants are the exponential forgetting factor uy, the initial covariance 
matrix P(0), and the initial parameter estimates, {a,(0)} and (8. (0)}; 
the order of the predictive model n, the scaling factor Bos and the 
assumed system time delay k. 

A simple numerical example due to Astrom and Wittenmark [1] 
was selected in order to study the effects of using different design 
options and initial parameter values. The example is a second order, 


Single-input single-output stochastic system with a time delay of one 
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21 
sampling interval and a noise variable of N(0,1). The system equation 
is 

y(t) - 1.9 y(t-1) + 0.9 y(t-2) 


= u(t-2) - u(t-3) + e(t) - 0.5 e(t-1) (2.30) 


The corresponding pulse transfer function model is 


A(q7') y(t) = B(q7!) u(t-2) + c(q7!) e(t) (2.31) 
where 
A(q7') = 1 - 1.9 q7! +.0.9 q72 
Bq) = 1-q" 
] 


cq!) =1-0.5q 


Note that the polynomial B has a zero on the unit circle thus 
violating one of the assumptions in section 2.2.1. Note also that 


polynomial A contains a zero on the unit circle. 


2.5.1 Derivation of Minimum Variance Regulator 
The predictive model corresponding to the example of 


equation (2.30) is 
y(t#2) + a, y(t) + a5 y(t-1) = a fu(t) + B, u(t-1) + 8, u(t-2)] 


+ ¢(tt2) (2.32) 


The identity in equation (2.5) is 


Rye0.5a Peettisl oqo + 094°) Ziq} fq es Bg) (2.33) 
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Solving equation (2.33) gives 


Z(q7!) =] + 1.4 qu! 


G(q7') = 1.76 - 1.26 q™ 


and substituting k, and polynomials B, Z and G into the predictive model 


of equation (2.20) gives 
y(t+2) - 1.76 y(t) + 1.26 y(t-1) = u(t) + 0.4 u(t-1) 


- 1.4 u(t-2) + ¢(t+2) (2.34) 


Comparing coefficients in equations (2.32) and (2.34) gives 


ay = -1.76 Boos 1.0 
aD = 1.26 | By = 0.4 
Bo = -1.4 


The minimum variance cuncrui i1aw for the system of equation (2.34) is 


u(t) = -1.76 y(t) + 1.26 y(t-1) - 0.4 u(t-1) + 1.4 u(t-2) (2.35) 


2.5.2 Typical Simulation Results 


A computer program has been written to simulate the 
response of the system in equation (2.30) when the self-tuning regulator 
is used. In all the simulation runs, the noise variable was specified 
to be N(0,1) using the IBM 1130 random number subroutine GAUSS which is 
also available on the IBM 1800 computer. The same random noise sequence 
was used in all runs. There were no constraints applied to the control 
Signal and R was set equal to the correct value of 1.0. For comparison 


purposes, the control law (2.28) was written as 


u(t) = & y(t) +... +8, y(t-m1) -8, u(t-1) - ... -B, (t-2) (2.36) 
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where 


Typical simulation results are shown in Table 2.1 and figures 2.2 to 
al 

Run number 1184 shown in figure 2.2 was used as a base case 
for purposes of comparison. For the base case, an exponential forgetting 
factor of y= 1.0 was chosen so that there was no exponential weighting 
of data. The initial parameter estimates (at t=0) and input/output 
data for t < 0 were assumed to be zero. The covariance matrix was 
arbitrarily specified as the identity matrix. The correct scaling 
factors By = 1.0, was used and the correct system order and time delay 
were used as the assumed values. Figure 2.2 illustrates that the 
convergence properties of the parameter estimates are very good but 
that the final estimates are biased as expected. (The correct values 
of the parameters are given in Table 2.1 and the recorded parameter 
estimates at the end of all runsare listed both in Table 2.1 and above 
the estimate curves in figures 2.2 - 2.6.) 

When exponential weighting of past data in the least squares 
estimation is introduced, the parameter estimates fluctuate more as 
shown in figure 2.3 and the rate of convergence is similar to that in 
figure 2.2. In general the proper selection of u can be made by 
examining the parameter estimate plots since too much fluctuation in 
the estimates implies too small a value of y is used, and slowly con- 
verging estimates means a smaller value of y should be used. In the 


Simulation study it was found that decreasing u(from 1.0) resulted in 
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increasing values of the accumulated loss. Therefore from the stand- 
point of control, it would be better to use a value of 1.0 for up in 
this example. 

The parameter estimates converged to the true values for 
only a small range of initial estimates, {a,(0)} and {8.(0)}. For 
example, with initial parameter estimates of 2.0, the estimates con- 
verged to incorrect values which contained large biases. Figure 2.4 
Shows that for initial values of estimates not too far from the true 
values, the convergence properties are very similar to those of the 
base case. It seems that the best choice of the starting values is 
zero if there is no prior information about the actual parameter 
value. 

Different values of the initial covariance matrix P(0) were 
tried (1 < P(0) < 1000 I) but there were no significant differences 
in the final parameter estimates. However, large initial values resulted 
in large initial fluctuations in the estimates as can be seen in figure 
2.5. The best results were obtained with P(0) = 2. ine effect oF 
different assumed values of By was also investigated. It was found that 
parameter estimates would converge for values of By which satisfied 


B 
0.6 < — Stes. HOM ae 2.0, the estimates tended to converge slower 
] 1 


initially than in the other runs, as shown in figure 2.6. In general, 
as BS increases, the convergence of the parameter estimates becomes 


slower. 


A comparison of the accumulated loss functions for the 
minimum variance and self-tuning regulator is shown in figure 2.7. The 


figure indicates that the major difference between the two curves occurs 
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during the transient period when the parameter estimates are adjusting. 
After the initial period, the responses are very similar. The corres- 
ponding output variables shown in figure 2.8 also support this observation. 
Figure 2.8 also includes the output of the uncontrolled system. This 
type of drifting response can be explained by considering the uncontrolled 
system transfer function: 
1 - 0.5q7! 
y(t) = —-—— —__ e(t) (2.37) 
Sige 208) 41-05 <) 

Since the denominator contains a root on the unit circle, this results 
in the integrating type of response shown in figure 2.8. 

The noise variable and the two control variables for the self- 
tuning and minimum variance regulators are shown in figure 2.9. After 
the initial transient, the behaviour of the control variables for the 


self-tuning and minimum variance regulators are very similar. 
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CHAPTER THREE 
STOCHASTIC EVAPORATOR MODELS AND DIGITAL SIMULATION STUDY 


3.1 Introduction 

In this chapter, several stochastic pulse transfer function 
models are derived starting from the fifth-order state-space evaporator 
model used by Hamilton [28]. The different models result from the 
various assumptions that can be made concerning the noise variables. 
The actual evaporator system and the computer program used in this 
Simulation study are also described. 

The latter part of this chapter presents the simulated 
responses of the double effect evaporator system when it is subjected 
to process and/or measttrement noise and unmeasured step disturbances. 
In each case, the open loop response of the system is presented 
followed by the closed loop response when a conventional multi-loop 
control system is used. Finally the conclusions of the open and 


closed loop studies are presented. 


3.2 Description of the Evaporator System 


A simplified schematic flow diagram of the double effect 
evaporator is shown in figure 3.1. The reader is referred to the 
Nomenclature section for definitions of the symbols. The evaporator 
is a forward feed unit which operates at a nominal feed rate of 
5.0 lb/min of 3.2 percent aqueous triethylene glycol and is heated 
by 2.0 1lb/min of 250°F saturated steam. The first effect is a 


natural circulation calandria type unit with 32 tubes, each of which 
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is eighteen inches long and has a 3/4 inch 0:D. The second effect 
is an externally forced circulation, long tube unit with three 
6 feet long, 1] inch 0.D. tubes. The second effect operates under 
vacuum and utilizes the first effect overhead vapour to concentrate 
the first effect bottom stream. The final product is about 1.6 1b/min 
of 10 percent aqueous triethylene glycol at the normal operating 
conditions. 

The IBM 1800 digital control computer in the Department of 
Chemical Engineering is interfaced to the double effect evaporator 
through appropriate converters and transducers. During normal 
operation the evaporator is controlled by means of six single control 
loops and four cascaded control loops. This conventional multi-loop 
control strategy is implemented through the standard Direct Digital 
Control (DDC) package. The most important controlled variable is 
product concentration C2 but the liquid levels, W1 and W2, must also 
be controlled within operating levels. In the multi-loop control 
system, levels WI and W2 are controlled by manipulating the two 
bottom flow rates, Bl and B2, respectively and the product concen- 
tration C2 is controlled by adjusting the steam flow rate S into the 
first effect. Further details concerning the evaporator equipment, 
instrumentation, and control loops can be found in the theses of 


Newell [29] and Jacobson [30]. 


3.2.1 Stochastic State Space Model 


Several models of the double effect evaporator have 


been developed in previous investigations by Newell [29] and Wilson [35]. 
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They range from tenth-order non-linear state-space models to first- 


order transfer function models. 


used is the fifth-order, linear, stochastic state-space model 


developed by Hamilton [28]. 


In the simulation study, the model 


This theoretical model was derived from 


linearized material and energy balances and adequately describes the 


process. This discrete, stochastic model is of the form: 
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the 


Au(t) + ad(t) + rw(t) 
t) + v(t) 
state vector (5 x 1) 


control vector (3 x 1) 
disturbance vector (3 x 1) 
output vector (3 x 1) 
measurement noise vector (3 x 1) 
process noise vector (6 x 1) 


sampling instant 


and 9, Lee and H are constant coefficient matrices of the 


appropriate dimensions. 


The state, control, disturbance and output vectors are 


defined as: 
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F! Wi! 
dite CF! Uae W2' 
HF C2! 


where W1', S', F' etc. denote normalized perturbation variables, i.e. 


Jn Fes (ss denotes the nominal 
: steady state value) 
Ss 


The process variables are defined in the Nomenclature section and 

their normal steady state values and the numerical values of the 
coefficient matrices are presented in the Appendix of Chapter Three. 

The measurement and process noise vectors v and w are assumed 

to be zero mean, uncorrelated, Gaussian noise sequences of N(0, OK) 
and N(0, Ob)» respectively. The six process noise variables in the 
vector w are the noise components of the three disturbance variables 
and the three control variables. Thus, coefficient matrix I is 


]. 


specified as I = Lm, 


He 


3.2.2 Derivation of the Stochastic Pulse Transfer Function Models 
Since the minimum variance and self-tuning regulators 
are based on stochastic models of the form of equation (2.1), such 
models will now be derived for the double effect evaporator. The 
starting point in the derivation is the stochastic state space model 
in equations (3.1) and (3.2). 
In order to derive a transfer function model of the evaporator, 
one output, one control variable and one noise element have to be 


selected from the fifth-order state-space model. Then the transfer 
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functions relating these variables must be derived in order to 
determine polynomials A(q@'), B(q7!) and C(q7!). 
The pulse transfer function matrix relating the output and 
control variables of the state-space model is given by 
H dj[I -1 -] 
-},_ Hadj[1 - sq Jaq 
Gdetes (atc andar ina (3.3) 
det | I= oq © | 
where “adj" denotes the adjoint matrix and "det" denotes the determinant. 
The pulse transfer function matrix relating the output and process noise 
variables is given by 
; -145..-1 
H adj[I - oq  Irq 


-] 
ge SS (3.4) 
=yw det | I - oq! | 


Once the output, control, and process noise variables are 
selected, then the corresponding pulse transfer functions are the 
corresponding elements of the matrices in equations (3.3) and (3.4). 


Let these pulse transfer functions be denoted by 


5 ba + boa +... 4 bea” q7'B(q7!) 
0 ee Oy (3.5) 
yu -] -5 -] 
T+ajq +... + aq A(q ) 
and 
af dq" + doa 5 Aes dq” D(q7!) 
Qy(Q) = eq oi) 
y Le ayq + eee + Gq A(q_) 


Thus when measurement noise is not present (i.e. v = 0), the output 


equation will be 


y(t) = dyy(a7 u(t) + aya”! )w(t) (3.7) 
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or 
-] -1 
y(t) = BOQ) yey + MG) wee) (3.8) 
A(q ') A(q ) 


where y, u and w are the selected elements of y, u and w, respectively. 
Note that equation (3.8) is in the form of the Astrém-Bohl in 
model in equations (2.1) and (2.3) with k=0 but d, = 0, and c, in 
equation (2.1) is equal to one. However, this difficulty can be 
easily resolved as shown below. 
Since {w(t)} is a sequence of independent N(0, on) random 
variables, then it follows that the delayed sequence, q7 w(t), has 


the same properties and can be denoted by e(t); i.e. 
q7lw(t) = e(t) (3.9) 


Let D(q7!) in equation (3.6) be rearranged as 


| d d 
D(q7!) = da7'(1 + ao toe ae ao) (3.10) 
Then equation (3.8) becomes 
-| d d d 
(t) = B(q_) (t-1) + u (1 + Ae fa} 2 q7*)q7 w(t) 
Ra Weer. carey) oi," q 
eae 
or ; 
- d d d 
y(t) = ee u(t-1) + (1 + = 4 Sa Tees = q”")e(t) 
A(q ) A(q ) 1 1 mie 


Comparing equations (3.12) and (2.3), it can be seen that a 
single-input, single-output stochastic model of the evaporator with 


only process noise has been derived with 
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= 


< Upon t 
A = dy C2 = q, 
wees d 
wa fer) Pr C, = d 
| 4 dy 
Cc, = 3 
2 qd, Cp = 0 


Alternative Models 

Two alternative transfer function models can be derived for 
the evaporator depending on the type of noise elements which are 
retained from the state-space model. 

In the above derivation e(t) was specified to be an element 
of w(t). However, if instead, e(t) was chosen to be the element of 
the measurement noise vector corresponding to the chosen output, 
then A(q7!) and B(q7!) remain unchanged but c(q7!) = 1]. This model 


would correspond to the predictive model that was used in section 2.4. 


A third type of model can be derived by retaining the 
appropriate element of v(t) in equation (3.2) during the derivation 


of the pulse transfer functions. This would give 


2 -1 
espe Peal abe Baty we s5y(t) (3.13) 
q 


where v(t) is the element of v(t) corresponding to the choice of y(t). 
If the additional assumption is made that v(t) = w(t) = e(t), then 


equation (3.13) becomes 


-] -] 
y(t) = Bia) u(t-1) + Mig —) e(t) (3.14) 
A(q ) A(q ) 
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m(q7!) = a(q7!) + d(q7!) (3.15) 
or 


1+ ma” ert dies mq” (3.16) 


a 
— 
2 
— 

it] 


Equation (3.14) can be written as 


= 
y(t) = Bigee). u(t-1) + aT (Vat mar ee ace meq” e(t) 
: (3.17) 


Again comparing equations (3.17) and (2.3), it can be seen that 


another Astrom-Bohlin model of the evaporator has been derived where 


A = 1 
c. =m, Pa ie” 
Note that this last model could also be derived by cisuming that 


the original state space model is of the form 


x(t+1) = ox(t) + au(t) + nd(t) + Ev(t) (3.18) 


y(t) = Hx(t) + v(t) (3.19) 


where r is a (5 x 3) matrix. That is, by assuming that the process 
and measurement noise vectors are identical. 

For the evaporator system, the most important control loop 
is the steam/product concentration loop, and the process noise 
variable of interest is the one corresponding to the feed flowrate 
since this is normally the most severe disturbance. Therefore a 
pulse transfer function model involving these three variables will 


now be derived. 
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The pulse transfer function relating product concentration 
C2 and steam S was obtained from the state-space model using the 
GEMSCOPE computer program [31] and the state-space model in the Appendix 


for Chapter Three. Since in the state space model, C2 is Y3> and S is 


Uy> this transfer function will be denoted by Ivy where: 
31 


(3.20) 


The pulse transfer function relating product concentration 


C2 and Wy» the noise term corresponding to feed flow F is 


(are) 


Therefore if measurement noise is not present, the stochastic pulse 


transfer function model for the evaporator is 


= -2 
3 pee BY UU Cd A ee 
y(t) 0.014 - 0.002 0.009 5 u(t-1) 


1-2.320q7! + 1.710q7* - 0.388q7 


eat -2 
0.002 + 0.0002 - 0.001 e(t) (3.22) 


1-2,320q7! + 1.710q7* - 0.388q72 
Note that all three polynomials A(q7!), B(q!) and c(q7!) in 
equations (3.20) and (3.21) contain a twice repeated zero of q=l. 


These factors cancel out and do not appear in the model of equation 
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3.2.3 Minimum Variance Regulator for the Evaporator 


For the evaporator model in the previous section, the 


time delay is zero. Thus equation (2.6) with k=0 reduces to 
-] 
AC MBE 


and from the identity in equation (2.5), 


] 2 


e(q7!) = 2.427 - 2.293q7! + 0.388q7 


Then the minimum variance controller for the model (3.22) is 
2.427 - 2.293q7! + 0.388q7> 
ut) = = 7 =, y(t) (3.23) 
0.014(1-0.127q7' - 0.679q~*) 


which can be written as 
ult) = ~-177s2ey{t).+ 167.4 y(t-1) - 28.3 y(t-2) 


+ 0.127 u(t-1) + 0.679 u(t-2) (3.24) 


3.3 Computer Program for the Simulation Study 

The objective of this simulation study was to investigate 
the response of the stochastic state-space evaporator model for open 
loop conditions and for cases when conventional multi-loop controllers 
were used. A computer simulation program called STR-2 was written 
for use on the IBM 1800 computer in the Department of Chemical 
Engineering. As shown in the flow chart of figure 3.2, the program 
begins with constant initialization and various option specifications. 
The outputs are calculated from the stochastic fifth-order state-space 
evaporator model in equations (3.1) and (3.2). Process and/or measure- 
ment noise of different noise levels can be added. Control options 


include no control, conventional multi-loop control, the self-tuning 
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Figure 3.2 Flow chart for computer program STR-2 
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regulator, or minimum variance control. For all closed loop runs, 
the control signals are constrained to conform to the physical 
limitations of the actual evaporator system. 

In the simulation study, evaporator variables of interest 
and parameter estimates of the predictive model for each run were 
printed. The evaporator variables and the parameter estimates were 
also punched on cards for the purpose of plotting. Two plotting 
program coreloads, RBNO2 and JCH10, written by Newell [32] and 
modified slightly by Hamilton [33] were used to plot the evaporator 


variables, 


3.4 Open Loop Results 

The open loop runs were made with six process noise variables 
each at the 10% level (i.e. Opn = 0.1) and with or without unmeasured 
step disturbances. These runs were made to give some indication of 
the fluctuations of evaporator variables. These runs were also made 
to serve as a comparison for the closed-loop runs. Different Gaussian 
random noise sequences with zero means were generated for the six 
process noise variables in each run and these same noise sequences 
were repeated for all other runs. The disturbance and manipulated 
variables that are plotted do not include process noise (except for 
the first figure), for the sake of clarity. This avoids the extensive 
overlapping of noisy curves that would result if the noise components 
were included. In all the simulation runs the system was assumed to 
be at the normal steady state at the start of a run. A list of the 


runs presented in this chapter is shown in Table 3.1. 
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TABLE 3.1 


OPEN LOOP AND MULTI-LOOP PROPORTIONAL CONTROL RUNS 


Note: 

All runs included six process noise variables (zero mean and 
standard deviation of 0.1) 

Run number 3259 included measurement noise in all three outputs 


(zero mean and standard deviation of 0.05) 


STEP ‘ 
FIGURE RUN NO. DISTURBANCES Keo ACCUM. LOSS (at t=300) 
er ee Ms ee ee 
3:3 3107 a - 0.0325 
3.4 3141 (-20%F, +20%F) a4 0.3805 } OPEN LOOP 
3.5 3255 (-302CF) : 6.1353 
36 3137 - 4,89 0.0325 
Py, 3136 é 0.0 0.0942 
3.8 3215 20.0 0.0217 Hae Eat 
: . ieee : PROPORTIONAL 
3.9 3130 (-20%F, +20%F)  -4.89 0.1215 | contRoL 
3.10 3176 (-302CF) -4.89 0.1061 


3.11 3259 io -4.89 0.2971 
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3.4.1 Effect of Process Noise 
As shown in figure 3.3, when the only disturbances are 

the six process noise variables with py = 0.1, all of the state 
variables remained fairly constant except for Wl, the level in the 
first effect. This shows that the magnitude of the process noise 
was not large enough to cause significant fluctuations in the evaporator 
variables. 

The manipulated and disturbance variables (S, Bl, B2, F, CF, TF) 
which are Gaussian random noise sequences are only shown in figure 
3.3. In subsequent simulation runs, these additive noise components 


are omitted to avoid overlapping of curves. 


3.4.2 Effect of Unmeasured Step Disturbances 


Figures 3.4 and 3.5 illustrate the effects of unmeasured 
step disturbances on the open loop system. In figure 3.4, a 20% step 
down in feed flow followed by a 20% step up approximately 50 minutes 
later were applied. While most of the variables remain within reasonable 
limits, the first effect level was drained completely due to the inte- 
grating nature of the liquid level. (The legend for these and succeeding 
figures is explained in the Appendix for Chapter Three. ) 

In figure 3.5, a 30% step down in feed concentration was 
applied. As a result the concentration in the first effect dropped 
which in turn caused the product concentration to decrease. 

Thus for the uncontrolled evaporator system, unmeasured step 


disturbances produce poor control. 
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3.5 Conventional Multi-Loop Control 


A simple control strategy that has been implemented on the 
actual evaporator system with fairly good success is multi-loop 
digital proportional control. As mentioned in section 3.2, the three 
major control loops are the two level loops and the product concen- 
tration loop. The primary control objective is to maintain product 
concentration at the normal steady state value. 

In this section, simulation runs were made using conventional 
proportional controllers in the three major control loops. The purpose 
of these runs was to investigate the effects of process noise, process 
and measurement noise and also the effect of unmeasured step distur- 
bances on the conventional control system. The nominal proportional 
gains were specified to be the constants used by Oliver [33]. The 
effect of using different proportional gains for the product concen- 


tration controller was also investigated. 


3.5.1 Magnitude of Proportional Gain in Steam-Product 
Concentration Loop 


The nominal proportional gains are [34], 


W1/B1 loop : Ky = 3.52 
W2/B2 loop : Kye = 15.8 
Ceise 1000p) 4: Koo = -4,89 


The control signals from the three loops were constrained to conform 
to the limitations of the real system. They were constrained as 


follows: 
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Figure 3.6 illustrates the response of the evaporator when 
multi-loop proportional control is used and six process noise variables 
are present. The control of C2 is fair while W1 is significantly 
improved over the open loop response in figure 3.3. 

Figures 3.7 and 3.8 show the evaporator responses when Koo is 
0 and -20.0, respectively. The levels, W] and W2, and the two bottom 
flow rates, Bl and B2, are omitted since they are very similar to the 
curves in figure 3.6. From figures 3.6 to 3.8, it appears that 
increasing the proportional gain Keo results in higher frequency 
oscillations and smaller deviations in C2. However when Kao is 
increased to -50.0 (figure not shown), the frequency of the C2 
oscillations increases but the magnitude of the deviations from steady 
state also increases. Thus the C2 response can be improved by 
increasing the numerical values of the gain only up to a certain limit, 


beyond which the performance deteriorates. 


3.5.2 Effect of Step Disturbances 


When a disturbance of 20% step down in feed flow 
followed by 20% step up in feed flow was applied, the C2 response 
in figure 3.9 initially increased and S decreased because of the 
sudden decrease in feed flow rate. This behaviour is characteristic 


of proportional control since there is always an offset after a step 
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disturbance. Furthermore, since a linear model is used in the study, 
the effects of step disturbances and process noise are additive. This 
is apparent from a comparison of figures 3.6 and 3.9. 

Figure 3.10 shows the evaporator response to a 30% step 
down in feed concentration. Once again the process noise variables 
give a certain pattern to the C2 curve and the step disturbance 


results in a general decrease in C2 and an offset. 


3.5.3 Effect of Measurement Noise 
When measurement noise with Omn = 0.05 was added to 

each of the three outputs and process noise is also present, the 
results in figure 3.11 were quite similar to the corresponding run 
without the measurement noise (see figure 3.6) except there were more 
fluctuations in the outputs. The measurement noise also produced 
large fluctuations in the three control variables S, Bl, and B2 as 
would be expected. Note that the output variables plotted in figure 
3.11 do not include the additive measurement noise elements, for the 
sake of clarity. 

Several other runs were tried using different proportional 
gains in the S/C2 loop. As the gain was decreased the C2 response 
improved. In fact, the best C2 response was obtained using zero 


gain in the S/C2 loop. 
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3.6 Conclusions 

Several stochastic pulse transfer function models of the 
evaporator have been derived starting from a stochastic state-space 
model. The particular form of the pulse transfer function depends 
on how the noise variable is specified (i.e. process noise vs measure- 
ment noise). The open loop simulation results indicate that 10% 
process noise in the three disturbance and three control variables 
does not significantly effect the evaporator. However, unmeasured 
step disturbances do have detrimental effects on the system. For 
example, a 20% step down in feed flow causes the first effect to 
drain completely. 

When conventional multi-loop proportional control is used, 
the process noise caused small fluctuations in all three outputs. 
With the addition of a step disturbance, larger offsets occurred. 
When measurement noise is added, the fluctuations in the actual 
product concentration increase while the three control signals 
exhibited much larger fluctuations. Increasing the proportional 
gain of the product concentration controller results in increased 
oscillations but smaller amplitude fluctuations in product concen- 
tration. However, this is only true up to a certain point, beyond 


which the control of C2 deteriorates . 
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CHAPTER FOUR 
A SELF-TUNING REGULATOR FOR THE DOUBLE EFFECT EVAPORATOR: 
SIMULATION RESULTS 


4.1 Introduction 

This chapter presents simulation results for the evaporator 
when the self-tuning regulator (STR) is used to control product 
concentration C2. It was mentioned in the previous chapter that it 
would be more realistic to use the fifth-order state-space model 
in the simulation rather than the pulse transfer function model 
because the actual evaporator is an interacting multivariable system. 
Therefore it would not be realistic to consider only a single-output, 
a single-input, and are roieo warjable and neglect their interactions 
with the other evaporator variables. 

For purpose of comparison with the conventional multi-loop 
control scheme, the two level control loops were left unchanged but 
the proportional controller for the S/C2 control loop was replaced 
by the self-tuning regulator. Thus the input variable of the pre- 
dictive model is steam flow rate, the output variable is the product 
concentration, and the noise variable is in effect a combination of 
all process noise and measurement noise variables plus the disturbances 
due to interactions between the evaporator variables. A block diagram 
showing the implementation of a self-tuning regulator in the multi- 
variable evaporator system is shown in figure 4.1. 

This chapter investigates the various design options that 


are available in the self-tuning regulator and its operation under 
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Figure 4.1. Block diagram of the closed loop system: the self-tuning 
regulator is used to control C2 and conventional 
proportional controllers are used to control W1 and W2 
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realistic conditions. Practical considerations in applying the self- 
tuning regulator and the conclusions of the simulation study are 
presented at the end of the chapter. A list of simulation runs for 


sections 4.2, 4.3 and 4.4 is shown in Table 4.1. 


4,2 Base Case Conditions 

The self-tuning algorithm has been explained in detail in 
section 2.4. Basically the algorithm can be divided into two parts 
at every sampling instant: estimation and control. Before the self- 
tuning algorithm can be used, the constants and initial values listed 


below must be specified: 


Exponential forgetting factor y 
Initial covariance matrix  P(0) 
Initial parameter estimates (0) 
Scaling factor Bo 

Model order n 


Time delay k 


The simulation study of the second order numerical example 
in section 2.5 has provided experience in the selection of these 
factors. Based on this experience a set of constants and initial 
parameters was selected for the base case run and served as a basis 
of comparison for the other runs. 

The exponential forgetting factor » was chosen to be 1.0 
so that there would be no exponential weighting of past data. The 
initial parameter estimates, {6,(0)}, were chosen to be zero based 


on the assumption that prior information about the evaporator system 
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was unknown. An initial covariance matrix of P(0) = 10,000 I was used 
since there was little confidence in the initial parameter estimates. 
Results obtained from the derivation of the pulse transfer function 
model were used to select the other constants for the base case run. 
The correct order of the predictive model and the correct time delay 
were used, that is, n=3 and k=0. Also the correct scaling factor, 
es 0.014, was used. In the simulation runs, the noise covariance 
was assumed to be the correct value, i.e. R = ObN This choice was 
made since preliminary runs indicated that the value of R was not 
critical. For example, R=] and R=0.01 gave almost identical results 
for the case when OPN = 0.01. Wittenmark and Wieslander [11] have 
recommended using R=1 if the actual noise covariance is not known. 
The base case run was carried out using these constants 
and initial parameters. The response of the system and the parameter 
estimates when process noise is present are shown in figures 4.2 and 
4.4, respectively. The large "bump" at the beginning of the C2 curve 
can be ecniained follows. Since zero initial parameter estimates 
were used, it takes several sampling intervals before the recursive 
least squares estimates have numerical values other than zero, and 
several more intervals before the estimates are near their final 
values (see figure 4.4). Therefore during this initial transient 
period, the control of C2 may not be good, which is the case in this 
run. The algorithm was able to adapt fairly rapidly, and after about 
20 minutes, the control of C2 was quite good. The control variable 
S exhibits bang-bang behaviour due to the large parameter values 


used in the control law. In figure 4.4 and succeeding plots of 
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parameter estimates the final values of the parameter estimates are 
noted in the figures. 

The parameter estimates in figure 4.4 of the base case run 
show large adjustments during the transient period. This was due 
to the large diagonal elements in the initial covariance matrix P(0). 
The rate of convergence of the parameter estimates was quite fast. 
Since the evaporator variables were generated from a fifth-order 
state-space model with six process noise variables, therefore, not 
only are the parameter estimates biased, they are also be 
different from the transfer function parameters derived in section 
3.2.2. Thus it is impossible to compare the estimated parameters 
from different runs with the derived parameters. Note also that 
different operating conditions or use of different initial parameters 
in the self-tuning regulator resulted in different parameter estimates 
at the end of each of the runs. 

Several different sets of process noise sequences were tried 
for the same base case run to check the randomness of the noise 
sequences. The evaporator responses (figures not shown) were very 
similar except for the initial "bump" of the C2 curves. Different 
process noise sequences appear to effect the size of the initial 
"bump" of the C2 curve and the one shown in figure 4.2 is the worst 
one. The final parameter estimates for different noise sequences 
were very similar. 

It was explained in section 2.4 that it may not be possible 
to estimate all the model parameters under closed loop conditions. 


To avoid this difficulty, By is usually assumed to be a known constant 
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as suggested by Astrom and Wittenmark [1]. Figure 4.3 shows a run 
in which all of the parameters including 8, were estimated. The 
evaporator response of this run is very similar to the base case 
run except the C2 curve contains no initial "bump" and has small 
sustained oscillations. The parameter estimates in figure 4.5 show 
that the {a} are very close to those obtained in the base case run 
and Bo is very close to the derived value in the pulse transfer 
function evaporator model in section 3.2.2. However, By> and Bo 
converge to very small values. Therefore it appears that even 
though the control of C2 is not very good, Bo can be estimated 
fairly accurately. Thus for some systems, on-line estimation during 
closed-loop conditions is one way of getting an initial value of By 
that can then be held constant for later runs. When successful, 


this approach eliminates the need for a priori specification of Bo: 


4.3 Initial Constants and Parameters for the Self-Tuning Algorithm 
The purpose of the simulation runs in this section is to 
study the response of the system when different initial values are 
specified for the self-tuning algorithm. The range of initial values 
which can be used without unduely upsetting the system or resulting 
in unstable behaviour was also determined. This is quite important 
because unless the self-tuning regulator can be used for a relatively 
wide range of constants and parameters, it will not be that attractive 
as a control strategy for use in industrial processes. 
In general, the selection of u, P(0), and {8.(0)} is not 


very difficult. Based on his simulation and experimental experience, 
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Wittenmark [7] has reported some guidelines on their selection. The 
value of up should be chosen to be between 0.99 and 1.0, and unless there 
is prior information available, 6,(0) = 0 should be used. If 6.(0) = 0 
is used, this meanslarge elements of P(0) should be used because of 

the low confidence in 6. (0). A range (101 < P(0) < 10,000 I) is 
recommended by Wittenmark [7]. These three parameters can be further 
adjusted by examining the response of the controlled system and the 
parameter estimate curves. The effect of different choices of these 
parameters are examined in sections 4.3.1 and 4.3.2. 

The selection of the constants for the regulator is more 
difficult since this depends on the user's knowledge of the physical 
system to be controlled. Therefcre, before the self-tuning regulator 
can be used, a rough idea of the system order and the process time 
delay should be known. As for Bo? according to Wittenmark [7] a 
wider range of values can be used if the control signal is constrained 
as was the case in the evaporator simulation. In sections 4.3.3 to 
4.3.5,the performance of the self-tuning regulator is evaluated for 
conditions where the system order, time delay and values of Bo are 


only approximately known. 


4.3.1 Exponential Forgetting Factor u 


The evaporator responses remained relatively unchanged 
for values of up in the range, 1.0 >.u 2 0.9, but the parameter 
estimates in figures 4.6 and 4.7 become more oscillatory as » decreases. 
When yp is decreased to 0.8, (figure not shown) the parameter estimates 


fluctuated widely and the evaporator response showed signs of becoming 


n , 
4 oe ? 
Vi 
‘ ce f. 
== i 
Bae 
i 
= * ~~ | 
ir rr Ail 
x W = 
"y ¢ r 
F ay < * es 
’ = > 
os — 
44 ve « 
x me. t23c ne e 
; i 4 
\ A a as 
i ieee - 
~ eae Se, 
‘i e « 
: 
“3 “at 
* 32 ’ ° * 
~ - 
3 x & 
a = , 
‘ey 
ri ~ ¢ 
* = . 
— 
™~ 
“ = ra , = 
Pee e 
va y ~ 
‘: : s 
Pate “ + . 
2 r eg => > £> . 
~ he 
a 
ie 1 
i oh a ae 1 2 dhe | 
2 re e177 Foo 6 
F , 
-~ ~~ 
; 7 & 5f . ltr | 
“YY Ga ae ae | 
Ais s y 


a: 
7“ 


_ ee x ce A ra _ a ‘a a 
‘dh . wa oe 2 7 Cal 
00 GAY tF boéuced, gine genie tee 
me 


7 at Bob w 


06°0 = 1 


*(YLS/WIS) SepPeULZsA uazeweued L*y7 aanbe4 96°0 = % S(YLS/WIS) SezPeWLISA UsazoWeUed Q9*p aunbly 
SSLANIN NI SAIL é S3LMNIN NI BAIL 
fe 8 u 8 N, 2 
8 8 8 g 8 8 


2 2 
¥ 3 
2 - 
a} Gj} 
5 : 5 
260°0 = “9 o L00°0-= “9 a 
: poet ii 
99/0 = 'Y ee p88'0 = *8 a 
SSINNIN NI SAIL SALPNIN NI SAIL 
f 8 >) 8 Mo ° i 8 x 8 u ° 
8 g 8 3 g g 2 g 8 3 8 g 
i 000-S- 
(O52 
bo 6] ee | 
= mo 
4A 8] 
5 O2L°0 .=~9 oe & 
in a 


™ 
Lo 
foe) 
— 

1 

" 

—_ 
3 


wt 
oO 
N 
oO 
" 
N 
3 
an SES 


v c00’S 


ae we 5 on 


ry ee 


- 
‘ ad ; 
r 4 * oe 
‘ | . 
. =a 
‘s . ony : 
‘ 
iga~e pt ; 

~\s Fe 4 3 ie : _ 
| | wea hal 

» \ y 7 | 

4 P 
- « : v . : . 
| | a ' Hf, ihe 

| * 5 ¥y 

I Peat, 
‘ ‘ 
; “ 
4 f 
‘ 
: 


-~ 


69 


unstable. In general, it is advisable to use a value slightly less 
than one for up [7]. This would ensure better tracking of model 
parameters especially parameters that are slowly time varying. 
Therefore the ability to track parameters and the quality of steady 


state control would depend on the user's choice of u. 


4.3.2 Initial Covariance Matrix P(0) and Parameter Estimates 
(6, (0)} 

The choice of the initial covariance matrix P(0) 
depends on the user's confidence in the initial parameter estimates 
{6,(0)}. For 6. (0) = 0, several values of P(0) were tried 
(10,000 I 2 P(0) 2 0.001 1). Typical evaporator responses are 
shown in figures 4.8 to 4.11. Reducing the elements of P(0) 
results in smaller elements in the gain vector K(t) which in turn 
causes smaller changes (per step) in the parameter estimates (see 
equations 2.25 to 2.27). This is quite evident when the two 
estimates plots are compared with the base case. For P(0) = 1001, 
the "bump" at the beginning of the C2 curve in figure 4.8 was reduced 
due to the smaller diagonal elements, but these elements were still 
large enough to ensure convergence of the estimates. For P(0) “1, 
the parameter estimates converged much slower and at the end of the 
run, the gain vector K(t) was too small to cause further changes in 
the estimates. Consequently, the control of C2 was not as good. 

The performance of the self-tuning regulator when non-zero 
initial parameter estimates and P(0) = 10,000 I are used is shown in 


figures 4.12 to 4.15. For both 6. (0) 25 and 6, (0) = =10, the 
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convergence properties were very good due to the large value of P(0). 
The C2 responses of the two runs were different during the initial 
transient due to the different signs and magnitudes of the initial 
parameters estimates. After the initial period, the control of C2 
was quite good in both cases. (Note that the scale of the oF plot 
in figures 4.14 and 4.15 was increased by a factor of two in com- 


parison with previous figures.) 


4.3.3 Scaling Factor By 


A large range of By values varying from 0.00001 to 1.0 
were considered. For 0.1 > Bite 0.00001, the product concentration 
was maintained within acceptable limits although the control of C2 
deteriorated at the two extremes, 0.1 and 0.00001, and there appears 
to be an optimum value of By in between. For Bo values outside 
this range, large sustained oscillations in C2 occur. When 
values close to Bai 0.014 were tried, the C2 responses were good 
and were very similar except for the initial "bump". Though a 
detailed search for the optimum value of B, was not performed, the 
best run, in terms of accumulated loss was observed when By = 0.03. 
Again, the lower accumulated loss is primarily due to the smaller 
initial "bump". 

In general, for relatively large values of Bo» the CF curves 
had large fluctuations and were still adjusting at the end of the run 
but the Bs curves converged quickly. The opposite situation occurred 
for relatively small values of 8): the oF curves converged but the Bs 


curves fluctuated. Typical results are shown in figures 4.16 to 4.19. 
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Wittenmark [7] has suggested that the user can examine the response of 
the control variable to find out if a good choice of By has been 

made. If the control signal is hitting the constrained limits too 
often, this means that By is too small; if the control signal is not 
hitting the constrained limits at all, this means that By should be 
decreased. These situations were observed by Wittenmark [7] in his 
simulation study and are also illustrated by the steam response in 
figures 4.16 and 4.17. Different values of Bo were also tried for 
smaller elements in P(0) and similar results were obtained. 

: If no control constraints are applied, a value of Pe 0.02 
results in unstable responses. This supports Wittenmark's comment [7] 
concerning the effect of control constraints on the operable range of 
Bo: 

4.3.4 Model Order n 

In the base case run, the correct model order of n=3 
was specified (see figure 4.2). Figures 4.20 and 4.21 show evaporator 
responses when n=2 and n=4 were used as the assumed model orders. The 
order of the predictive model did not seem to influence the C2 response 
very much after the initial transient period. Note that the "bump" in 
the initial part of the C2 response increases as the assumed model 
order increases. This can be explained as follows. When a higher 
Rete ddel is used, more parameters have to be estimated and since 
zero initial conditions were assumed, the information vector y» contains 
only zero elements initially. If the recursive least squares algorithm 


in section 2.4 is examined closely, it can be seen that the elements of 
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the parameter vector have non-zero values only if the corresponding 
information vector also has non-zero elements. Therefore, for higher 
order models, it takes more sampling intervals before the information 
vector and the parameter vector have non-zero elements. As a result, 
higher order models will generally give poorer control during the 
initial transient period. 

Additional runs were performed with n=1 and n=5 (figures not 
shown). For n=], the regulator has only one parameter and is actually 
a proportional controller; the evaporator response was similar to 
those shown in section 3.5. For n=5, the evaporator response was 
Similar to the run where n=4. The parameter estimates for the runs 
where 1 < n < 5 all had good convergence properties. Figures 4.22 


and 4.23 show the parameter estimates plots for n=2 and n=4, respectively. 


4.3.5 System Time Delay k 


Simulation runs ne also made using assumed time delays 
of k=] and k=2. Since the actual time delay in the derived pulse 
transfer function model is k=0, control of C2. became worse as k 
increased, as would be expected. However, for k=2 the control was 
still satisfactory which suggests that overestimating the delay may 
not degrade the performance of this system too much. Figures 4.24 
and 4.25 show the evaporator response and the parameter estimates for 
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4.4 Performance of the Self-Tuning Regulator Under Non-Ideal 
Operating Conditions 


An evaluation of the performance of the self-tuning regulator 
for a variety of operating conditions was a prime objective of the 
simulation study. In this section the effects of process noise level 
and the inclusion of measurement noise and unmeasured step disturbances 
are described. A list of simulated runs for this section is shown in 


Table 4.2, 


4.4.1 Effect of Process Noise Level, Opn 


In the base case run, a 10% process noise level was 
used (i.e. Opn = 0.1). Figure 4.26 and 4.27 show evaporator responses 
when 1% and 50% process noise levels are used, respectively. Increasing 
the noise level resulted in an increase in output variance and large 
fluctuations in all the other evaporator variables, as would be 
expected. Note that the parameter estimates in figure 4.28 and 4.29 
are quite similar although increasing the noise level results in more 


fluctuations. 


4.4.2 Effect of Measurement Noise 
Addition of 5% measurement noise to the three evaporator 
outputs resulted in poor control of C2 as shown in figure 4.30. The 
corresponding parameter estimates are shown in figure 4.32. A run with 
Only measurement noise and no process noise (not shown) gave results 
similar to figures 4.30 and 4.32. These results indicate that 5% 


measurement noise is more significant than 10% process noise. 
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It is interesting to note that the latter part of the C2 
response in figure 4.30 is worse than the response in figure 3.1] 


where proportional control is used with K -4.89 for the same 


CPs 
process conditions. 


4.4.3 Effect of Unmeasured Step Disturbances 


‘In practical applications it is seldom feasible to 
measure all possible disturbance variables. Consequently, the per- 
formance of the self-tuning regulator when unmeasured, non-zero mean 
disturbances occur is quite important. 

For a low process noise level of Opn = 0.01 and two step 
disturbances in feed flow, the self-tuning regulator results in 
satisfactory control as shown in figure 4.31. The corresponding 
parameter estimates in figure 4.33 show that the estimates change 
rapidly after the first step disturbance, but not after the second 
step change. This somewhat unexpected behaviour can be explained as 
follows. After the first step down in feed flow, the product concen- 
tration increases slightly (this can not be observed in figure 4.33 
but can be clearly seen in the computer printout). Therefore, the 
output elements in the information vector y (t) are above their steady 
state values, and this causes the elements in the covariance matrix 
P(t) to change. This adjustment in P(t) resultes in larger elements 
in the gain vector K(t) which in turn causes some adjustment of the 
parameter estimates. These changes in the estimates when used in the 
control law seem to have the same effect as increasing the gain in a 


proportional controller. This can be observed by comparing the steam 
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flow rate in figures 4.26 and 4.31. 

The set of parameter estimates at t = 75 min were able to 
handle the second step disturbance quite well. As a consequence, the 
product concentration, covariance matrix and gain vector were all 
observed to undergo much smaller changes which meant there was little 
change in the estimates. 

Figures 4.34 and 4.35 show the evaporator response for 10% 
process noise and step disturbances in feed flow and feed composition 
respectively. The C2 responses in both runs are very similar to the 
response for the base case conditions (see figure 4.2). These results 
indicate that the self-tuning regulator can handle unmeasured step 
disturbances quite well. The corresponding parameter estimates are 
shown in figures 4.36 and 4.37. Note that the estimates in figure 
4.36 do not undergo sudden changes after the first step change at 
t = 25 min, as was the case in figure 4.33. This was attributed to 
the higher process noise level in figure 4.36 since the 10% noise level 
tended to reduce the effect of the 20% step changes in feed flow rate. 
It is interesting to note that C2 in figure 4.35 does not have an offset. 

When all the parameters of the predictive model (including 8.) 
were estimated and unmeasured step disturbances were applied, the control 
of C2 was not very good. Figures 4.38 and 4.39 show two cases where 
feed flow and feed concentration disturbances occur. The two C2 
responses were significantly worse than the same case without the step 
disturbances (see figure 4.3). This illustrates that poor C2 responses, 
can result if all parameters are estimated during closed loop operation. 


The corresponding parameter estimates plots are shown in figures 4.40 
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and 4.41. For all runs where Bo is estimated, the (8,} converged to 


values very close to zero. 


4.5 Practical Strategy for Implementing Self-Tuning Regulator 


When a process control computer is used to control an 
industrial process, it will probably not be economical to replace 
most of the conventional digital controllers in individual loops 
with self-tuning regulators. This is because self-tuning regulators 
require more execution time due to the recursive least squares equations 
and they also require more oo race Space for input/output data. The 
process computer might not be able to handle the extra work load or 
have sufficient storage space for the extra data. One way of getting 
around this problem is to use the self-tuning regulator only on 
important control loops and/or after process conditions have experienced 
a significant change. Furthermore, the self-tuning regulator need 
only be applied until good parameter estimates and a satisfactory 
controller have been obtained. Then the estimation step can be by- 
passed and the now constant parameter estimates can be used in the 
control calculations (i.e. a controller with constant coefficients in 
the control law). 

Figure 4.42 shows a run in which the controller had constant 
parameters which were obtained from the final parameter estimates at the 
end of the base case run. Even when two 40% step disturbances in feed 
flow were applied, the controller was able to control C2 quite well. 

Theoretically, if the estimator is being used for very long 


periods of time the elements in both the covariance matrix and the gain 


~ - 7 
f : ie > 32434 2 
rs ¥ 
oe 
, % - i on @ > 
3. 365 PAT aS 
ale af Sora oma 
a ; - - 
i « le +) “— “ nV: “ur 
ee pee SAY Vary © tog 
r aa owe : 
: -~ > ha 
5 en . i , P. - 
a —— , 2 7 
Tt a y 7 a r fated hn 
Crs 7, tan = c SSG boRX ret wide 
z a ‘Ss as + 9 
= _ os - : 
: ‘ s o wemalyg re 
= ‘2 oo ee oe j ‘a AL 
b 5 7 +m ts <- 
_' ~ 
* 
“ i, ~~ 
+ od auh ae? col 
Ses cen - 
‘7 — —s 
: F 
* = 4 at ~ 
. Pg ~ att © 
= 
a 5 7 sag oye 5 
v ‘ ~ ay —— aft 


( A: of mol dova. 2 i 
i ’ J 


‘ 


y 
or 2% | 
\pa2e. 


98 


Joze,nBau ul suazaweued quezsuod *(YIS/4 %Ov+ 64 %OV-/WIS) Bsuodsau sozeuodergz 


Zp'by aunbe4 


SALNNIN NI SIL SALMNIN NI SAIL 


OOT SZ OS S2 0° Set OOT Sc 


y 


NIWG1 NI 401d WSLS 


NIW/G1 NI 3ivevO1d 


OS 


RZ 


bt fe 


OcH ‘NI NI “TAD 


Te 


4 v 


DIA 3d NI NIN 


Ot 


99 


vector in the algorithm decrease as the parameter estimates converge. 
This is satisfactory as long as process conditions do not change. But 
if the process conditions do change, the elements of the gain vector 
are now too small to be of any effective use in updating the parameter 
estimates. To avoid this type of situation, the regulator should be 
re-initialized periodically. This can be done by re-initializing the 
covariance matrix. Figure 4.43 shows a run in which the initial parameter 
estimates were the final estimates obtained from the base case run and 
P(0) = 100 I was employed at the beginning of the run. The corresponding 
parameter estimates are shown in figure 4.45. Note that the estimates 
deviate from the initial values at the beginning of the run but gradually 
approach these values towards the end of the run. 

One way of avoiding the problem of poor estimates and poor 
control during the transient period has been suggested by Wittenmark [7] 
and is as follows. Suppose a well-tuned conventional controller is used 
initially and least square estimation is performed to determine the 
parameters of the predictive model. Then the switch to the self-tuning 
regulator can be made after the parameter estimates have become fairly 
constant. 

Figures 4.44 and 4.46 show a run in which a proportional 
controller was used in the S/C2 joop for the first 14 minutes with 
Keo = -4,89. The self-tuning regulator was then substituted for 
the Proportional controller at the 15 minute mark. Note that the 
initial "bump" in previous C2 curves is eliminated since the self- 
tuning regulator starts with better initial parameter estimates, 


namely, those that were generated during proportional control. Thus 
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this strategy provides a systematic technique for a bumpless transfer 


from a conventional control strategy to the self-tuning regulator. 


4.6 Conclusions 

From the extensive simulation study of the double effect 
| evaporator model, it can be observed that the self-tuning regulator 
is very simple to use. The initial choice of constants and parameters 
required for the regulator can be made without detailed knowledge of 
the process to be controlled and the regulator performs fairly well 
even for poor choices. Once the regulator is operational, these 
constants can be adjusted to better values by examining the process 
data and the parameter estimates. 

The self-tuning regulator also performed well when different 
unmeasured disturbances of relatively large magnitude were applied, 
but did not fare as well when 5% measurement noise was added. 

Practical ways of implementing the self-tuning regulator 
were suggested in section 4.5. These methods tend to ensure good 
performance of the regulator when it is used periodically on a long 
term basis. 

In general, the self-tuning regulator definitely performed 
better than the conventional proportional controller in the simulation 
study. This success results because the algorithm is designed to 
minimize the output variance in the presence of noise. 

The simulation results showed the good convergence properties 
of the self-tuning algorithm. The system was stable unless very poor 


constants were used in the algorithm. The algorithm itself is very 
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simple and there are no complex mathematical operations involved. 
The compactness and the simplicity of the algorithm plus its ability 
to handle relatively large process noise levels and unmeasured 


disturbances make it very attractive for industrial applications. 
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CHAPTER FIVE 
SELF-TUNING REGULATOR FOR THE DOUBLE EFFECT EVAPORATOR: 
EXPERIMENTAL RESULTS 

5.1 Introduction 

This chapter presents experimental results of the application 
of conventional multi-loop control and the self-tuning regulator to the 
double effect evaporator. The pilot plant evaporator that is interfaced 
to the IBM 1800 process control computer has been described in Chapter 
Three. The simulation results for the evaporator have been presented in 
Chapters Three and Four. This chapter also includes a brief description 
of the control computer program that was used. The experimental results 
are then presented and compared with the simulation results. Finally, 
conclusions meine practical uses of the self-tuning regulator in 


industry are given. 


5.2 Experimental Procedure 


In order to compare the simulated and experimental results, a 
Similar set of runs was performed on the actual evaporator system. 
Initial conditions were maintained close to those used in the simulation 
runs. Prior to the start of each experimental run, the evaporator was 
brought to the nominal steady state conditions by means of multi-loop 
digital proportional-integral] controllers (the nominal steady state 
values are shown in the Appendix for Chapter Three). The control program 
was then initiated and the evaporator data were collected via 
individual DDC (Direct Digital Control) loops at every sampling interyal 
and were stored on disk files. The control and sampling periods were 


specified to be 64 seconds as in the simulation runs. The control 
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Signals calculated from the control algorithms were set as setpoints to 
the individual DDC control loops. At the end of each run, all the 
collected data were printed and also punched on computer cards. Documen- 
tation of a typical experimentation run is shown in the Appendix for 
Chapter Five, including heat and material balances and a log of pressures, 
temperatures and flow rates etc. These records plus the computer print- 
out of the evaporator variables and parameter estimates, and the 
recorder charts for each run have been compiled in an “Experimental Run 
Data Book". 

A statistical analysis of normal steady state experimental 
data has been performed by Hamilton [28] to estimate the actual process 
and measurement noise levels. These estimates were found to be much 
less than the 10% and 5% levels used for process noise and measurement 
noise, respectively, in the simulation study. The actual estimates 
calculated by Hamilton were 0.1% to 5% for the process noise and 0.1% to 
2% for the measurement noise. In all the runs, the evaporator variables 
were filtered by simple filters such as exponential or Union filters so 
that smooth and accurate data would be used and recorded. The three 
output variables were filtered by exponential filters with filter 


constants of 0.75. 


5.3 Computer Programs for the Experimental Study 


The real-time computer programs that were used in the closed- 
loop experimental study were basically the ones developed by Newel1 [36]. 
Program coreloads RBN30, RBN31, RBN32 and RBN33 were modified so that 
three different control strategies could be used: multi-loop propor- 


tional control, the self-tuning regulator and minimum variance control. 
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To avoid confusion with the original program coreloads, these modified 
coreloads were re-named, FLC30 to FLC33. The initialization procedures 
of these control programs are essentially the same as those used by 
Newel] and these procedures can be found in Newell's "User's Manual" [36]. 
Programs RBN1I9 and RBN2O, used to print and punch the accumulated 
experimental data, were also modified to include the parameter estimates 
and were re-named as FLC19 and FLC20. Plotting programs RBNO] and RBNO2 


[32] were used to plot the evaporator variables. 


5.4 Experimental Evaporator Results 


To verify the main results of the simulation study, a similar 
set of experimental runs were performed. Again, conventional multi-loop 
control techniques such as digital proportional control and proportional 
plus integral control were applied. Then the self-tuning regulator was 
ig in the steam/product concentration control loop and the effects of 
various design decisions and parameter values were inyestigated, as in 
the simulation study. 

In addition to the lower process noise levels in the actual 
evaporator system, several other changes had to be made for the experi- 
mental runs. First, tighter constraints were placed on the steam flow 
rate which was limited to between 0.5 and 2.5 Ib/min. The upper bound 
was chosen to limit the pressure in the first effect to a safe operating 
level. The lower bound was chosen to be 0.5 1b/min instead of zero to 
prevent negative pressure in the first effect as a result of stoppage 
of boiling and condensation of overhead vapour. Second, the proportional 
gains in the two level control loops had to be reduced to less than 


half the normal values which were used in the simulation study. The 
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EXPERIMENTAL RESULTS FOR CONVENTIONAL CONTROL STRATEGIES 


TYPE OF 
FIGURES RUN NO. CONTROL 
Sel. es 26 ML prop. control 
Bee 54 2/7 ML prop. control 
5.5 21 ML PI control 
v0 24 ML PI control 
CET | 25 ML PI control 


STEP 
DISTURBANCES 
(420% F, -20% F) 
(-30% CF) 
(+20% F, -20% F) 
(+20% F) 
(-20% F) 
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reason for reducing the gains was that the bang-bang characteristics of 
the steam signal resulted in large fluctuations in the first effect 
pressure and the first effect level measurements, which in turn affected 
the first effect bottom flow rate, and the second effect level and bottom 
flow rate. These interactions between process variables resulted in 

poor control of product concentration. To avoid this undesirable 
situation, the proportional gains in the two level control loops were 
decreased in order to reduce the interactions. 

These differences in operating conditions should be kept in 
mind when the simulation results are compared with the experimental 
results. The experimental runs for conyentional multi-loop control and 
the self-tuning regulator are shown in Tables 5.1 and 5.4, respectively. 
Since the experimental runs were performed for different periods of 
time, the parameter estimates listed in Table 5.4 and in the figures 


were those recorded at the end of the runs. 


5.4.1 Conventional Multi-Loop Control 


Conventional multi-loop proportional control using the 
modified version of the multivariable control programs [36] was applied 


with the following proportional gains for the three control loops. 


TABLE 5.2 
PROPORTIONAL CONTROLLER GAINS 


S/C2 loop ;: Koo = -4,89 
W1/B1 loop 5 Ky ae 
W2/B2 loop . A Ke = 6.5 
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Note that the proportional gains in the two level control loops are lower 
than those used in the simulated study of Section 3.5. 

The evaporator responses for feed flow and feed composition 
step disturbances are shown in Figures 5.1 and 5.2. The C2 responses 
and the two effect levels and bottom flow rates in both figures are less 
oscillatory then those in the simulation results. This is due to the 
lower process noise level and lower gains in the level control loops. 
Both the feed and first effect concentrations (CF and Cl) were only 
measured before, but not during, the experimental runs, thus, the Cl 
curves are not shown and CF is a straight line in Figure 5.1. 

During these closed loop runs, the least squares estimator in 
the self-tuning algorithm was used to estimate the parameters of the 
predictive model. The initial constants and parameters specified for 
the algorithm were those used in the base case run in Section 4.2. The 
parameter estimates for these two runs are shown in Figures 5.3 and 5.4. 
The convergence of the estimates is not very good because the evaporator 
variables undergo little excitation. 

The DDC version of multi-loop proportional plus integra] 
control was also applied since this is a widely used conyentional 
control technique. The closed-loop response of the evaporator to two 
20% step disturbances in feed flow for King and McNeill‘s DDC 
controller constants [37] in Table 5.3 is shown in Figure 5.5. (KP and 
KI are the proportional and integral constants, respectively, readers 
are referred to the DDC Manual [39] for more information on DDC control 


loops). 
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TABLE 5.3 
PROPORTIONAL AND INTEGRAL CONTROLLER CONSTANTS 


CONTROL KING & MCNEILL [37] HAMILTON [38] 


LOOP KP KI KP KI CONTROL INTERVAL 

S/C2 o20 0.04 Salen 0.0234 64 sec 
W1/B1 6.0 0.009 6.0 0.009 2 sec 
W2/B2 32.0 0.009 32.0 0.009 2 sec 


The control of C2 was better than in Figure 5.1 as expected, and the 
level control was also improved. The DDC controller constants used by 
Hamilton [38] shown in Table 5.3 were also tried and the resulting 
evaporator responses are shown in Figures 5.6 and 5.7. Comparing 
Figures 5.5 and 5.6, it can be seen that King and McNeil1's constants 
resulted in smaller fluctuations in Bl and B2 but larger fluctuations in 
S. Also the control of C2 after the first step disturbance was better. 
From these runs, it appears that conventional multi-loop 
strategies provide fairly good control of C2 and also keep the two levels 


near their steady state values even when step disturbances are applied. 


5.4.2 Self-Tuning Regulator 


The experimental runs for the self-tuning regulator are 
summarized in Table 5.4. 
The base case run in Section 4.2 was also tried experimentally 
and the response is shown in Figure 5.8. Process yariables C2, WI, W2, 
Bl, B2 and S all have sustained oscillations even though no step 


disturbance was applied. This can be explained as follows. At about 
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the 10 minute mark, C2 decreases because of the lower steam flow. As a 
result, the steam flowrate increases and boils off more water from the 
first effect thus lowering WI. Because of the decrease in WI, Bl 
decreases (due to the W1/B] control loop) which in turn causes W2 and B2 
to decrease. (C2 increases due to S and Bl, and the self-tuning regulator 
then reduces the steam flow rate to compensate for this. When the steam 
Signal drops, the whole procedure is repeated but in the opposite 
direction, From that point onwards, all the key variables have sustained 
oscillations. The corresponding parameter estimates plots are shown in 
Figure 5.10 and all have good convergence properties. Again, the final 
parameter estimates are listed in the figure. 

Several unsuccessful attempts were made to eliminate these 
oscillations. Values of By of 0.03, 0.05, and 0.005 were tried but the 
results were similar to those in Figure 5.8. Then the control and 
sampling periods were reduced by half but the oscillations still 
occurred. The two level proportional controllers were then replaced by 
the standard PI controllers (with Hamilton's constants), but this did 
not provide any improvement as can be seen from Figure 5.9. The 
corresponding parameter estimates are shown in Figure 5.11. 

By examining the results obtained when self-tuning regulator 
was used, it seemed that the strong interaction initiated by the bang- 
bang behaviour of steam could be reduced by reducing the proportional 
gains in the level control loops. Thus, a run was tried using zero 
gains in the level controllers. The result was that Bl and B2 were held 
constant and the control. of C2 was good (figure not shown), but the two 
levels dropped below their operating ranges due to the integrating nature 


of the two holdups. To avoid this situation, another run was tried 
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using very small gains in the two level control loops. The self-tuning 
regulator was able to keep the evaporator at the nominal steady state 
as shown in Figure 5.12. When 20% step disturbances in feed flow were 
applied, the self-tuning regulator was able to handle the disturbances 
without difficulty as shown in Figure 5.13. Note that with the low 
gains used in the level loops, the increase in feed flow resulted in a 
large upsurge in the first effect level. The corresponding parameter 
estimate plots for these two runs are shown in Figures 5.14 and 5.15. 
The convergence of the estimates was quite good. 

In order to prevent large offsets in W1] after a step disturb- 
ance in feed flow, the proportional gains in the level loops were 
increased to the values shown in Table 5.4. These gains which are 
about half the values used in the simulation study were also used in the 
remaining experimental runs. Two runs were performed using the base 
case parameters in the self-tuning algorithm with two different types of 
step disturbances as shown in Figures 5.16 and 5.17. In Figure 5.16 the 
C2 response was good and the regulator was able to handle the two 20% 
step disturbances in feed flow quite well. With the larger gains in the 
level control loops, the levels were controlled much better. The C2 
performance was comparable to the analogous run in the simulation study 
(see Figure 4.34). In comparison with the conventional P and PI control- 
lers (see fig. 5.1 and 5.5) the self-tuning regulator produces smaller 
deviations in C2 but resulted in small sustained oscillations towards 
the end of the run. In Figure 5.17, the 30% step disturbance in feed 
concentration had little offset on C2 as compared to the proportional 
control run in Figure 5.2. The strong influence of S on the levels and 


the bottom flow rates is again apparent in Figures 5.17. The correspond- 
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ing parameter estimates are shown in Figures 5.18 and 5.19. 

Figures 5.20 and 5.21 illustrate the effect of using values of 
By = 0.1 and By * 0.001 in the self-tuning algorithm. The characteristics 
of the steam signals are similar to those observed in the simulation 
study (see Figures 4.18 and 4.19). For the larger scaling factor of 
By = 0.1, the steam signal exhibits similar small fluctuations initially 
but it eventually results in sustained oscillations. For the smaller 
scaling factor of Pn 0.001, the steam signals hit the constraints 
more frequently. The parameter estimates in Figures 5.22 and 5.23 show 
even a closer resemblance to those obtained in the simulation study. 

For ia 0) .eche fa} have large fluctuations and the {83} converge, 
while for sei 0.001, the exact opposite can be observed. 

Model orders of n = 2 and n = 4 were used in the self-tuning 
regulator and gave the results shown in Figures 5.24 and 5.25, respect- 
ively. Step disturbances in feed flow were applied to both runs and 
the control of C2 was good. For n= 2, the C2 performance was one of 
the best observed in the experimental study. The corresponding 
parameter estimates are shown in Figures 5.26 and 5.27, respectively. 

Experimental runs using the self-tuning regulator in which all 
the parameters including 8, were estimated, proved to be unsuccessful 
as compared to those runs in the simulation study. The steam signal 
would drop to the lower limit and stayed there until the runs were 
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5.4.3 Practical Strategies for Implementing the Self-Tuning 
Regulator 
Practical strategies for applications of the self-tuning 
regulator were considered in the simulation study of Section 4.5. In 
this section, two of these strategies were tried experimentally. 

In Figure 5.28, the parameter estimates recorded at the end of 
the experimental run in Figure 5.14 were used as constant parameters in 
the control law of Equation (2.28). This constant parameter control law 
is the minimum variance controller for the empirically determined 
transfer function model. The satisfactory C2 response verifies that 
once good parameter estimates are obtained, the estimator can be by 
passed and these constant parameters can be used in the control law 
until conditions warrant the use of the estimator again. 

The strategy of using the estimator to estimate model 
parameters (except By) while the evaporator was under conventional 
proportional control was also tried. The results are shown in Figure 
5.29. The estimator was used for the first 24 minutes and the resulting 
parameter estimates were then used as the initial parameter estimates 
at t = 25 minutes when the conventional controller was replaced by the 
self-tuning regulator. A scaling factor of 8, = 0.02 and an exponential 
forgetting factor of u = 0.995 were used in this run. 

The purpose of this strategy was to generate good parameter 
estimates before the self-tuning regulator is initiated. In the 
simulation run it was found that this method eliminated poor control 
during the transient period when the parameter estimates were adjusting 
(see Figure 4.44). In all the experimental runs, because of the lower 


process noise level, the control of C2 was satisfactory during the 
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transient period, so the merit of using this strategy is not as 
apparent. The corresponding parameter estimates are shown in Figure 5.31. 
A similar run was performed using the same large level loop 
gains that were used in the simulation study and the base case initial 
parameters in the self-tuning regulator. The results are shown in 
Figure 5.30. The C2 response was good for the major portion of the run, 
but towards the end, the high gains in the level loops again caused 
sustained oscillations in the control and output variables. The 


corresponding parameter estimates are shown in Figure 5.32. 


5.5 Conclusions 

For multi-loop proportional control, the experimental responses 
contained less fluctuations than those obtained in the simulation study. 
This was due to the lower process noise leyels in the actual] evaporator. 
The experimental results also show that the steady state contro] of C2 
was fair and that C2 responses were satisfactory when step disturbances 
were applied. The standard DDC version of multi-loop proportional plus 
integral control gave slightly better control of C2, as expected, and 
the C2 response would probably be improved if there were looser control 
over the two effect levels (i.e., lower gains). All the conventional 
multi-loop control strategies that were applied gave much better C2 
responses than those reported by Jacobson [304s 

The experimental application of the self-tuning regulator for 
the base case conditions of the simulation study was unsuccessful. All 
of the control and output variables had sustained oscillations due to 
control loop interactions. However, when the gains in the leyel control] 


loops were reduced by a factor of two the self-tuning regulator applica- 
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tion was successful. With the smaller proportional gains in the level 
control loop, good experimental results were obtained for base case 
conditions. A direct comparison of simulation and experimental runs 
could not be made, since the proportional gains of the level loops were 
different, the steam signals were constrained within narrower bounds, 
the actual process noise levels were lower, and output variables contain 
small measurement noise. In general, the performance of the self-tuning 
regulator was quite good and the characteristics of the regulator 
observed in the simulation study were verified in the experimental 
results. The corresponding parameter estimates did not converge as well 
as in the simulation runs but were still satisfactory. 

Three different values of 8, were tried and the characteristic 
influence of B, over the steam signal and the parameter estimates was 
similar to the simulation results. Different model order were tried 
(n = 2,3 and 4) and all gave satisfactory C2 responses. These experi- 
mental results verified the observation in the simulation study that the 
choice of model order was not critical. 

It was also demonstrated that once a good set of parameter 
estimates have been obtained, these parameters could be held constant 
and used in the regulator until there were significant changes in 
process conditions. The proposed strategy of estimating the model 
parameters under conventional feedback control for a certain period 
before switching to the self-tuning regulator did not show any improve- 
ment for the actual evaporator system. This occurred because the 
evaporator was close to steady state during the estimation period, and 
with low process noise levels and little or no excitation in the steam 


signal, the estimator could not generate good estimates for use in the 
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self-tuning regulator. 

Experimental runs in which all the parameters (including By) 
were estimated were not successful in contrast to the successful 
simulation runs. In all of the unsuccessful runs, the steam signal 
would drop to the lower limit initially and remain there until the run 
was terminated. Thus it appears that attempting to estimate all the 
parameters in the self-tuning regulator results in non-identifiability 
of the parameters and poor control of C2. Consequently, one parameter 
should be fixed for the evaporator system as recommended by Astrém and 
Wittenmark [1]. 

To summarize, the self-tuning regulator performed quite wel] 
even with an approximate set of initial parameters and handled unmeasured 
step disturbances quite well. In comparison with well-tuned P or PI 
controllers, the self-tuning regulator resulted in similar steady state 
control but better control when unmeasured step disturbances occurred. 
However, the merits of the self-tuning regulator were not fully demon- 


strated in this application due to the low process noise levels. 
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CHAPTER SIX 
CONCLUSIONS 


The purpose of this investigation was to provide a detailed 
evaluation of the self-tuning regulator and to investigate its 
performance in a practical application to a pilot plant, double effect 

evaporator. Simulation and experimental studies have demonstrated that 
the evaporator application was successful and have shown that the 
regulator works well for a relatively wide range of design parameters. 
The following design parameters were investigated: exponential forget- 
ting factor, initial covariance matrix, initial parameter estimates, 
scaling factor, model order and system time delay. The self-tuning 
regulator was also evaluated for non-ideal operating conditions including 
uiirerent process noise and measurement noise levels, and unmeasured 


Step disturbances. 


The conclusions concerning the effects of the various design 


parameters can be summarized as follows: 


1. Scaling factor BOS The scaling factor is the most difficult para- 
meter to specify. In the simulation study a wide range of values 
(0.1 2 B, 0.00001) gave satisfactory results when the control 
signal was constrained. This was reduced significantly when the 
control signal was unconstrained. An experimental run with 

pars 0.1 resulted in large fluctuations in the output indicating 

that only a smaller range of values can be used in practice. Both 


the simulation and experimental results show that smaller values of 


By produce significant fluctuations in {Ba}; rapid convergence of 
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Car and large, frequent fluctuations in the control signal. The 
converse situation occurred for large values of Bo: Here {B.) 
converged quickly, {a.} fluctuated and the control signal was 


smoother. 


2. Model order n; The assumed order did not significantly affect the 
process responses. Except for the initial twenty minute transient, 
second, third and fourth order evaporator models gave practically the 
same simulated responses, but in the experimental study, under- 
estimation of the model order (n = 2) gave better control of the 


output. 


3. System time delay k; Overestimation of the process time delay by two 
sampling intervals (128 seconds) still resulted in satisfactory 
control in the simulation study. This factor was not investigated 


experimentally since the actual evaporator system has no time delay. 


4. Initial covariance matrix P(0); The selection of the initial 
covariance matrix and initial parameter estimates should be made 
together since the choice of the initial covariance matrix is based 
on the degree of confidence one has in the initial parameter estimates. 
If poor initial parameter estimates are used, then large diagonal 
elements in the initial covariance matrix should be assigned. Large 
djagonal elements in P(0) produce large initial adjustments of the 
parameter estimates followed by fast convergence. During the initial 
transient large variations in the control signal occur and unless the 
control signal is bounded, this will produce undesirable results. 


Smaller diagonal elements in P(0) produce a slower rate of convergence 
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and smaller changes in the parameter estimates during the initial part 
of the response. In the simulation study, choosing P(0) = 100 I 
provided good convergence of the parameter estimates starting from 
zero initial values. However, P(0) = 10,000 I was used in most simu- 


lation runs and in all experimental runs to ensure fast convergence of 


the estimates to satisfactory values. 


5. Initial parameter estimates 6(0); The self-tuning regulator worked 
well in the simulation study for a wide range of poor initial estimates 
including 6. (0) = 0, -5, and 10, when large diagonal elements were used 
in P(0). Zero initial parameter estimates were used in all experimental 
runs and rapid convergence of the parameter estimates occurred. Thus 
a priori knowledge of actual system parameters is not necessary for 


satisfactory performance of the self-tuning regulator. 


6. Exponential forgetting factor nu; Values of » between 0.9 and 1.0 
did not significantly affect the simulated responses but the 
parameter estimates fluctuated more for smaller values. A value 


of » = 1.0 proved to be satisfactory in the experimental study. 


In this application, it was found that the design parameters, 
particularly By? P(0), and u can be adjusted on-line to give better 
values after examining the transient process response and the parameter 
estimates. 

The simulation results showed that the self-tuning regulator 
can handle unmeasured step disturbances better than a digital propor- 
tional controller. The experimental runs confirmed this result but 


gave small oscillations around the steady state in some runs due to the 
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significant fluctuations in the manipulated variable. In the simulation 
study, the output variance increased as the process noise level 
increased and that the addition of 5% measurement noise resulted in 
large fluctuations in the controlled variable. It was found that 
estimating all of the parameters including the scaling factor was 
possible in the simulation study though the control of the output was 
not very good. However, this method did not work well in the experi- 
mental runs. 

The interacting nature of the multi-input, multi-output 
evaporator system strongly affected the performance of the self-tuning 
regulator. This was clearly seen in the experimental results where the 
influence of other evaporator variables, such as the bottom flow rates, 
resulted in large sustained c2cittations when normal proportional gains 
were used in the two level control loops. These interactions were 
significantly reduced by decreasing the gains in the level control loops 
but very small oscillations about the desired steady state still 
persisted in most runs. 

In general, the self-tuning regulator performed well especially 
when a good set of initial parameters were used. The regulator is very 
simple to use and does not involve complex mathematical operations. It 
can take into account system time delays and can have different numbers 
of control parameters depending on the assumed system model. It can be 
used continuously so that it will automatically adjust control parameters 
on-line to suit changing operating conditions or process dynamics, or it 
can be used periodically for the purpose of tuning control parameters. 


A practical scheme for implementing the self-tuning regulator in which 
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a conventional controller is used initially before switching to the 
regulator has proved to be useful in this work and is recommended. The 
self-tuning regulator is probably best suited for stochastic processes 
with relatively high process noise levels and frequent unmeasured distur- 
bances since its design is based on stochastic control theory. Its use 
is recommended especially if conventional control techniques are unable 


to provide good control for these conditions. 


6.1 Future Work 
The simulation and experimental studies in this thesis provide 

a detailed examination of the self-tuning regulator and general 

guidelines for its use. However, there are a number of areas that 

require further investigation and these will be discussed as follows. 

1. The self-tuning algorithm can be extended to include feedforward 
control of a measured disturbance [7,22]. This extension can be 
applied to the evaporator system by using measurements of either 
feed flow rate or feed concentration. Although it has been shown 
that the self-tuning regulator can handle unmeasured disturbances 
well, this feedforward inyestigation may prove valuable for other 
systems. 

2. The self-tuning regulator can also be modified to include setpoint 
control of the output variable [7]. This extension would increase 
the capability of the regulator. 

3. The scaling factor in the regulator was normally assumed to be 
constant in this study. The effect of holding some other parameter 
constant instead of the scaling factor or the effect of fixing 


several parameters could also be studied. 
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Only a single self-tuning regulator was used for one particular 
control loop in the evaporator application. A logical extension 
would be to use self-tuning regulators in several control] loops 
either sequentially or simultaneously. 

In this study, the self-tuning regulator was used in an interacting 
multivariable system and the effects of interaction were not 
directly taken into account. Howeyer, from previous evaporator 
studies, it is known that the first effect bottom flow rate affects 
the product concentration to a significant extent. The influence 
of this control variable could be included in the self-tuning 
regulator by using a two input, one output process model. This 
extension may reduce the experimental difficulties caused by control 
loop interactions. 

The self-tuning regulator is usually designed to minimize the 
output variance. Other performance indices that include control 
variable terms have been recently proposed [15,16,23] and should be 
evaluated. These modifications might reduce the relatively large 
magnitude fluctuations in the control variable. 

The extension of the self-tuning type of control strategy to 
multi-input, multi-output systems is currently under inyestigation 
[21]. If successful, such an extension would certainly increase 
the flexibility of the control strategy and the evaporator would be 


an appropriate system for an application. 
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NOMENCLATURE 


Polynomial in Astrém - Bohlin model 
Coefficient in polynomial A 
Polynomial in Astyom - Bohlin model 
Coefficient in polynomial B 

First effect bottom flowrate 
Second effect bottom flowrate 
Polynomial in Astrom - Bohlin model 
Coefficient in polynomial C 

First effect concentration 

Second effect concentration 

Feed concentration 

Polynomial 

Coefficient in polynomial D 
Disturbance vector 

Expectation operator 

Random normal variable 

Feed flow rate 

Polynomial 

Coefficient in polynomial G 

Output coefficient matrix 

First effect enthalpy 

Feed enthalpy 

Identity matrix 


Gain vector 
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N 
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Time delay 

Number of 8 parameters 
Polynomial 

Coefficient in polynomial M 
Normal distribution 

Order of model 

Overhead from first effect 
Overnead from second effect 
Covariance matrix 

Pressure in first effect 
Pressure in second effect 
Backward shift operator 
Covariance of noise yariable 
Steam flow rate 

Temperature in first effect 
Temperature in second effect 
Feed temperature 

Sampling instant 

Control vector 

Loss function 

Measurement noise vector 
First effect holdup 

Second effect holdup 

State yector 

Output yector 

Polynomial 


Coefficient in polynomial Z 


152 


a“ ed 
» ‘J rr 
ai lk ee ae 
tess ia ass ti WIOH 
¢ Vs 


Sra ep 
r i ‘ 
yo ee 


‘ = by ‘ 
rt 
x. a . 
NORS2 TTT PES Tay 
viatim sons itavod 
+ 
oT a eee ae 
‘ = Wt ~& 64% t 
— f F - 
. toe? o292 af ayers 
‘ : m? _ 
Tari nyt ‘tile, bumpatoe 0B 
"he ve eT 


oe 


* 
3 


ey , ‘ hes 
> Ss Poe 


Ch j Tan} 6 Ps r IRS ty Cv 


uh Bess: T 
4 aah Sa 
Pare auth ar OFimy 
- ls “i 
‘2% bagee at. sete 
set iok 
' yf ses \y ae 


= Mal 

? 

® 2a ts a a> 

,.* 
ai 


b) 


c) 


Greek 


Coefficient of predictive model 
Coefficient of predictive model 
Scaling factor 

Process noise coefficient matrix 
Control coefficient matrix 
Moving average process of driving noise e 
Parameter vector 

Constant 

Exponential forgetting factor 
Disturbance coefficient matrix 
Standard deviation 


Information vector 


Subscripts 


jth element 


Vector 

Matrix 

Process noise 
Measurement noise 


Steady state 


Superscripts 


Matrix transpose 


Matrix inverse 


Estimated value 


Perturbed variable 


e) Abbreviations 

adj Matrix adjoint 

det Determinant 

prop Proportional 

CC Concentration controller 

CR Concentration recorder 

FC Flow controller 

FR Flow recorder 

LC Level controller 

PC Pressure controller 

KP Proportional constant 

KI Integral constant 

PI Proportional plus integral 

f) Codes for Computer Plots 

Paes Denote time of step disturbance 

<jrs Denote initial steady state 
Type of Run ;: 

SIM Simulated 


EXP Experimenta] 
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Disturbance (+ / X% / XX): 


ee hag Positive or negative step 
X% Step size as percentage of steady state 
XX Process variable disturbed 


Control Mode: 


OL Open loop 
ML Multi-loop conventional control 
STR Self-tuning regulator 


PI Proportional plus integral 
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APPENDIX FOR CHAPTER THREE 


(a) Normal Steady State Operating Conditions 


x 


|< 


{oO 


Ro 


Five element state vector 


Wl First effect holdup 

Cl First effect ccneeniration 

Hl First effect solution enthalpy 
W2 Second effect holdup 


C2 Second effect concentration 


Three element control vector 


S Steam flowrate to first effect 
Bl First effect bottoms flowrate 


B2 Second effect bottoms flowrate 


Three element disturbance vector 


F Feed flow rate 
CF Feed concentration 


HF Feed enthalpy 


= [W1, W2, C2] 


16] 


46 1b (16.0 in) 
4.59% glycol 
190 BTU/1b 

42 1b (11.0 in) 
10.11% glycol 


2.00 1b/min 
3.49 1b/min 
1.58 1b/min 


5.0 1b/min 
3.2% glycol 
162 BTU/1b 
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(b) Matrices for the Discrete Evaporator Model (based on Steady State 


condition in (a) and a sampling interval of 64 seconds) 


-0.0008 -0.0912 0 0 
0.9223 0.0871 0 0 
-0.0042 0.4377 0 0 
-0.0009 -0.1052 1.0 0.0001 
0.0391 0.1048 0 0.9603 


-0.0119 -0.0817 0 
0.0116 0 0 
0.1569 0 0 

-0.0137 0.0847 -0.0406 
U.UIS/ -0.0432 0 

| & 


0.1182 0 -0.0050 
-0.0351 0.0785 0.0049 
-0.0135 -0.0002 0.0662 

0.0012 0 -0.0058 
-0.0019 0.0016 0.0058 
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APPENDIX FOR CHAPTER FIVE 
The following is a listing of computer printouts 
containing information about the controllers and steady state condi- 
tion of the evaporator that are collected at the start of a typical 
experimental run. The reader is referred to the User's Manual [38] 
for information on the computer programs used. 


(a) Computer printout from program FLC32 


eee SELF-TUNING CONTROL EXECUTIVE *x#% 


INITIALIZATION DATA FOR RUN 12 


DATA STORAGE INTERVAL(SEC)= 0 

COUNTER INDICATOR= 0 

CONTROL INTERVAL(SEC)= 64 

JFLAG= 2 KFLAG=99 

STD. DEV. (PROC. NOISE)= 0.10 STO. DEV.(MEAS. NOTSE)= 0.00 
KP(1)= -4.89 KP(2)= 1.50 KP(3)= 6.50 


INITIAL CONSTANTS 
BO=0.01400 M= 3 K=0O Z= 1.000 


INITIAL CONTROL PARAMETERS 
0.00 0.00 0.00 0.00 0.00 


INITIAL COVARTANCE MATRIX 


10000 .001 0.000 0.000 0.000 0.000 
0.000 19090 .001  .000 0.000 0.000 
0.0900 0.000 10000 .601 0.000 0.000 
0.000 0.000 0.000 10000.091 0.000 
0.000 0.000 0 .000 0.000 10000.001 


STEADY STATE VALUES FOR XM, Uy AND D 


O.16060E 02 0.22140E 03 OeLOS3Z0F OZ 0.99899E-01 
0.18009E Ol 0.31480E 01 0.13580E Ol 
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(b) Computer printout from program PVLOG 


THE FOLLOWING IS A LISTING OF THE MEASURED VALUE FOR EACH LOOP 
RECORN SPECIFIED, WITH THE ASSOCIATEN MEASUREMENT CONDITION CODES. 


15 APR 75 

15/27 HRS 
LOOP ID MEASUREMENT LOOP DESCRIPTION CONDITION CODES 

GTV ERR PVR) MEAS 
141 3.063 L8/MIN FEED WATER FLOW l 2 2 
139 1.433 LB/MIN FEED SOL. FLOW 1 2 2 
149 4.494 LB/MIN TOTAL FEED FLOW 1 2 iz 
151 3.143 LB/MIN Bl 1 2 2 
152 1.439 LB/MIN B2 i 2 2 
132 1.798 LB/MIN STEAM FLOW 1 2 2 
150 1.540 LB/MIN CONDENSATE FROM SECOND EFFECT i! 2 2 
126 1.385 LB/MIN CONDENSER CONDO. FLOW 1 2 2 
133 48.590 L8S/MIN- COOLING WATER FLOW 1 2 4 
136 16.070 IN H20 FIRST ESFECT, LEVEL 1 2 2 
137 11.900 IN 20 SEPARATOR LFVEL 1 2 2 
125 5.000 IN H20 CONDENSATE LEVEL 1 2 2 
135 4.810 PSIG FIRST EFFECT PRESSURE 1 2 2 
1324 -17.960 IN HG SECOND EFFECT PRESSURE a 2 2 
i142 190.300 DEG F FEED WATER TEMPERATURE 1 2 2 
138 190.000 DEG F FEED SOL. TEMPERATURE 1 2 2 
159 188.600 DEG F TOTAL FEED TEMPERATURE 1 2 2 
143 0.032 FEEN CONCENTRATION 1 1 1 
164 221.600 DEG F FIRST EFFECT TEMPERATURE ] 2 2 
129 0.099 ACTUAL PRODUCT CONC. l 2 2 


MEASUREMENT CODES USED ARE 


GTVLU ERROR CODES 


1 NO ERROR NETECTED 

2 VALUE NOT IN LOOP RECORD 

3 INVALID IT0. LOOP NOT IN TABLF 

4 INVALID DATA IN CALL GTVLU STATEMENT 


PROCESS VARIABLE RECORD (LOOP RECORD) CONDITION CODES 
1 LOOP OUT NF SERVICE 
2 LOOP IN SERVICE 


MEASUREMENT VALIJF CONDITION CODES 
l INPUT BAD 
2 INPUT NORMAL 
3 LOW ALARM LIMIT 
4 HIGH ALARM LIMIT 
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MPX 
~16318 
-16316 
-16314 
-1 6313 
=16312 
-16311 
-16306 
-16305 
-16304 
-16303 
~16302 
=1630 ft 
-16300 
lo. 09 
-16298 
-16297 
-16296 
=16295. 
~16294 
-16293 
-16292 
-16291 
-16290 
-16289 
-16288 
-16287 


-16286 
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Computer printout from program BJ20D 


STEADY STATE TEMPERATURES 


TEMP .(F) 
219.9 
218.9 
154.8 
153.0 
162.2 
219.2 

72.6 
160.9 
74.2 
13352 
190 .6 
189.6 
189.3 
85.5 
298.3 
253.1 
252.4 
221.4 
218.3 
188.5 
182.8 
eS 
Li2.9 
161.6 
161.7 
161.4 


194.0 


DESCRIPTION 
TR20 - FIRST EFFECT, SOLUTION AT TOP OF DOWNCOMER 
Thos = FIRSY EFFECT, SOLUTION ABNVE THE TUBES. 
TR14 — FIRST EFFECT, SOLUTION AT THE ROTTOM. 
TR33 - SECOND EFFECT, CONDENSATE AFTER THE COOLER 
TR 6 ~— SECOND EFFECT, SOLUTION AT THE BOTTOM, 
TR32 -— SECOND EFFECT, STEAM BLEED. 
TR27 - STEAM, OUTLET OF THE FEED WATER HEATER. 
TR34 - SECOND EFFECT, PRODUCT JUST BEFORE PUMP, 
TR21 —- SOLUTION TANK, TNLET SOLUTION. 
TR18 - 
WES SES, TOTAL AFTER MIXING. 
TT11 —- FEED, WATER BEFORE MIXING. 
TT19 - FEED, SOLUTION BEFORE MIXING. 
TT13 - CONDENSER, CONDENSATE INSIDE. 
TR15 — STEAM, TO FIRST EFFECT BEFORE ORIFICE. 
TR13 - FIRST EFFECT, STEAM IN THE CHEST. 
TR 5 - FIRST EFFECT, STEAM CONDENSATE. 
TR19 -— FIRST EFFECT, SMLUTION AT ROTTOM OF DOWNCOMER 
TR 2 -— FIRST EFFECT, VAPOUR LEAVING. 
TR 7 - FIRST EFFECT, FFED AT INLET. 
TR 4 -— SECOND EFFECT, FEED AT INLET. 
TR16 - SECOND EFFECT, STEAM AT TOP OF CHEST. 
TR 9 — SECOND EFFECT, SOLUTION AT TOP. 
TR17 — SEPARATOR, SOLUTION LEAVING. 
TR12 — SEPARATOR, VAPOUR LEAVING. 
TR10 - SECOND EFFECT, STEAM AT INLET TO CHEST. 
TR28 - SECOND EFFECT, STEAM CONDENSATE BEFORE COOLER 
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-16284 99.4 TR 1 — CONDENSER, COOLING WATER LEAVING. 
-16283 90.4 TR30 ~ RUN DOWN TANK, COOLING WATER LEAVING. 
-16282 99.1 TR31 — RUN DOWN TANK, LIQUID LEAVING. 

-16281 86.0 TR11 -— CONDENSER, CONDENSATE LEAVING. 


(d) Computer printout from program HMBAL 


PSnAPR 7S 
15/728 HRS 


ERROR OF CLOSURE ( PERCENTAGE) 
INPUT-OUTPUT 


OVERALL BALANCES 


TOTAL MASS BALANCE @.177 LB/MIN € 3-94 ) 
COMPONENT MASS BALANCE 2.223 LE/MIN ¢ cone, 2 
HEAT BALANCE @.122E 23 BTU/MIN ¢ 4.99 3} 


FIRST EFFECT BALANCES 


TOTAL MASS BALANCE -2.212 LB/MIN ¢ -4.66 ) 


SECOND EFFECT AND CONDENSER BALANCES 


TOTAL MASS BALANCE 2.387 LB/MIN ¢ 2.24 » 
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